Microscopic Theory for Multiple Light Scattering in Magnetic Fields 



Bart A. van Tiggelen 1 , Roger Maynard 1 and Theo M. Nieuwenhuizen 2 
1 Laboratoire de Physique Numerique des Systemes Complexes, 
Universite Joseph Fourier/CNRS, 
B.P. 166, 38042 Grenoble Cedex 09, France. 
2 Van der Waals-Zeeman Institute, 
Valckenierstraat 65 - 67, 1092 AN Amsterdam, the Netherlands. 

(June 27, 1995) 



Abstract 

Q\ We present a microscopic theory for multiple light scattering occurring in inhomo- 
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geneous 3D media subject to an external magnetic field. Magneto-optical effects (the 
Faraday effect and the Cotton-Mouton effect) occur inside the small inhomogeneities. 
We thereby take into account the spatial anisotropy, time-reversal-symmetry break- 

^ ^ ing and birefringence caused by the magnetic field, and discuss the consequences for 

CN the diffusion tensor and the polarization characteristics of the diffuse light. We will 

frequently compare our findings to a similar phenomenon in dilute polyatomic gases: 
^vq the Beenakker-Senftleben effect. Coherent Backscattering and the field-field corre- 

lator (Ei (B) Ej (0)) are addressed, which have both been obtained experimentally 
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by the group of Maret etal. All modifications in transport theory due to the mag 



netic field exhibit a rather sensitive dependence on the scattering phase shift of the 
individual scatterers. 
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I. PHYSICAL CONTEXT 



Multiple scattering of light is a fascinating topic having much overlap with other branches of 
physics such as solid-state physics, atomic physics and astrophysics. Most interesting aspects are due 
to the persisting role of interference in multiple scattering of light, together with the relative ease to 
observe and manipulate it in controlled laboratory experiments. For electrons this is more difficult due 
to phase-destroying mechanisms such as electron-phonon coupling that are not so easy to eliminate. 
The interference in multiple light scattering manifests itself in - by now rather well understood - 
phenomena as Coherent Backscattering [1] and universal fluctuations [2]. Other applications, such 
as localization of light, are still part of speculation. 

One perhaps unexpected conclusion is that interference in multiple light scattering is in fact very 
difficult to suppress. Experimentalists spend hard times in undoing the effect of speckles in light 
measurements and extracting an average value of some transport variable. Absorption is known to 
suppress multiple light scattering but not the relative importance of interference. A finite coher- 
ence length of incident light does not destroy all interference and some experiments, like Coherent 
Backscattering can even be done with sunlight. 

The only known way to manipulate interference externally is by applying a magnetic field. In the 
solid state a magnetic field breaks the charge symmetry between holes and electrons, destroying the 
Weak Localization correction and giving rise to a negative magneto-resistance. In the case of light, 
magneto-optical effects in the dielectric constant have a similar impact on Coherent Backscattering 
of light. The physics of the latter, constructive interference of time-reversed waves, is sometimes 
believed to be similar to the Weak Localization of electrons [3]. The dominant magneto-optical effect 
is the Faraday effect. Mathematically it is due to an anti-symmetric term in the dielectric tensor 
linear in the external field vector B [4]. In a homogeneous medium it is well known that this causes a 
rotation of the polarization vector of a linearly polarized plane wave proportional to traversed length 
and magnetic field. Equivalently, two opposite circularly polarized light beams achieve different 
speeds. As such it is often a mean to determine interstellar magnetic fields. 

The study of the Faraday effect in inhomogeneous media is much more recent. The inclusion of 
Faraday rotation into a transport equation of light has first been carried out by plasma astrophysi- 
cists [5] in order to understand the combination of emission and absorption of polarized synchrotron 
radiation in optically thick ionized plasmas. However, in that case the light scattering is inelastic. 
In sharp contast to the light processes in plasmas, elastic scattering in disordered dielectric media 
conserves the phase so that interference can take place. After initiating work by Golubenstev [6], 
MacKintosh & John [7] found, in their rather profound study, that the Coherent Backscattering 
of light will be suppressed by the time-reversal symmetry breaking character of the Faraday effect. 
As has been shown by Martinez and Maynard [8] [9], elastic multiple light scattering from inhomo- 
geneities in a dielectric background with Faraday active material is - under some approximations - 
extremely suited for numerical simulation. These simulations also work for the cases of higher mag- 
netic fields and the inhomogeneities replaced by genuine Mie spheres, both of them being difficult to 
attack analytically. 

Experimentally the Coherent Backscattering of light in magnetic fields is studied extensively by 
Erbacher, Lenke and Maret [10] [11]. Some of these results will be discussed from a theoretical 
point of view in section V. In these experiments it is difficult to separate between "medium" and 
"scatterers" . Nevertheless, roughly speaking, they can be divided into two kinds. In some of them 
the magneto-optical effects occur in the medium and not in the immersed particles. This is the 
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situation envisaged in the theory of MacKintosh & John and the numerical simulations of Martinez 
& Maynard. In other work, the magneto-optical effects occur inside the particles and not in the 
surrounding medium. For this case no theory exists and the present paper aims to fill this gap. We 
will show that this case allows - in the conventional Boltzmann approximation of transport theory - 
an analytical solution incorporating all aspects of the magneto-optical effects in multiple scattering 
of light. Most of them, like the Cotton-Mouton effect and the anisotropic diffusion tensor, have not 
been discussed before. 

In pursuing analogies with other branches of physics we wish to mention the Beenakker-Senftleben 
effect in dilute gases [12] [13] [14]. This effect embodies the field-dependent and anisotropic nature 
of transport coefficients for paramagnetic and even diamagnetic polyatomic gases in an external 
magnetic field. This is caused by the precession of the magnetic moment of the molecules between 
the collisions. This rotation tends to average out the nonspherical collision cross-section between the 
molecules, and thereby lowers transport properties like viscosity and heat conductivity. The magnetic 
influence on the scattering cross-section itself is known not to be responsible for the Beenakker- 
Senftleben effect, since it can be shown to be negligible [15]. 

In our case it is the electric polarization vector that rotates between the collisions of the light 
with the dielectric particles. If one assumes that the cross-section is unaffected by the magnetic 
field indeed a great resemblance with the Beenakker-Senftleben effect turns up. One property of 
this phenomenon, the presence of a one-parameter scaling variable u p t (uj p the precession frequency 
and t the mean free collision time), has also been seen in multiple light scattering experiments [11]. 
Although some differences will nevertheless show up, we will demonstrate that the basic transport 
variable for multiple light scattering - the diffusion constant - becomes anisotropic in a magnetic 
field, quite similarly as in the Beenakker-Senftleben effect. Even an anti-symmetric part of the 
diffusion constant appears, generating a transverse current, perpendicular to both density gradient 
and magnetic field [16]. An anisotropic diffusion tensor (signifying a transverse current) is best 
known for conduction electrons and ionized plasmas [5], where the Lorentz force gives rise to the 
Hall-conductivity (and a less well-known transverse heat conductivity called the Righi-Leduc effect). 
On the other hand, in previous treatments of magneto-optical effects, the direction of the magnetic 
field was not noticed to play a role in the diffusive regime. 

The outline of the paper is as follows. In the next section we derive the scattering properties of 
one magneto-optical particle and define most of our notation. In section III we obtain the average 
electric field amplitude and find the dispersion relation u (k, a) of a coherent beam with wave vector 
k and helicity a. 

In section IV we calculate the incoherent energy density by summing - following conventional 
techniques - the Ladder diagrams in the Boltzmann approximation. In the diffusion approximation we 
obtain the diffusion tensor for the light in the magneto-optical medium. The associated eigenfunction 
- the "Goldstone mode" - turns out to be anisotropic in polarization space. The translation of bulk 
results to finite media is crucial to compare to experiments. It will turn out that role of the "skin 
layers" [17] for incoming and outgoing light is very important. Due to the anomalous step length 
distribution in these skin layers, the Stokes parameters in diffuse transmission will not be equal to 
zero in a magnetic field. We will show that the magnetic field can induce oscillations in the Stokes 
parameters as a function of the slab thickness and the field strength. Similar oscillations have been 
reported on the base of numerical work [8], but have not yet been observed. 

In section V we will discuss Coherent Backscattering by considering the most-crossed diagrams. 
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We find the dephasing length due to the change of reciprocity relations (the Faraday effect). The 
absence of the usual reciprocity also causes the diffusion constant featuring in the most-crossed 
diagrams to be different from the one in the Ladder diagrams, as is also known to be true for 
electrons in the solid state [18] [19]. Attention will be paid to polarization properties of Coherent 
Backscattering, as well as to the anisotropy of the line shape for various directions of the magnetic 
field. Again, we shall rely on a strict diffusion approximation, and ignore (eight) non-diffusive 
eigenvalues of the most-crossed diagrams. So far, this has also been done in the interpretational 
analyses of the data [11]. Let us recall that the angular behavior in the center of the backscattering 
peak is due to long-range diffusion only. 

In section VI the correlation function (^Ei (B) Ej (0)^ is considered, that is the correlation func- 
tion of the electric field with and without external magnetic field (in transmission). In the C\ 
approximation this correlation function is related to the intensity correlator (7(0) 7(B)) which is 
experimentally obtainable [10]. We find a difference in phase between E (0) and E (B) (showing 
up as a complex eigenvalue near zero) giving rise to a phase factor in the field correlation function 
itself. The field correlation function {Ei (B) E 3 (0)^ has not yet been measured but recent advances 
in heterodyne detection methods may facilitate this. 

Although the model is perhaps too simplistic to provide a quantitative theoretical picture for 
present experiments, it is the simplest approach starting directly from Maxwell's equations. Along 
the way, no approximations other than the usual ones will be made (independent scattering approxi- 
mation, diffusion approximation). Therefore we hope that this paper will guide to interpret existing 
data qualitatively and perhaps will even initiate future experiments. 

We want to draw attention to Appendix A where we define our tensor notation. 



II. SINGLE SCATTERING 



In this section we derive the scattering properties of one magneto-optical dielectric scatterer 
located in vacuum. By taking the scatterer •pointlike we shall be able to obtain these scattering 
properties in analytic form so that they can be used for multiple scattering in principle without 
further approximations [20] [21]. In this model we are even able to add a scattering resonance, being 
characterized by a relatively large cross-section, together with a delay and out-of-phase response of 
the scattered wave. The phase function of this model is only influenced by anisotropy in polarization 
indices caused by the magnetic field, and not by anisotropy in wavenumbers as common for Mie 
scatterers. 

Throughout this paper we set c = 1/y/eo = 1. The refractive index - a tensor of rank two - of a 
particle in an external magnetic field is assumed to be given by, 

n (B) = ml+— <fr + MBB (1) 

to 

In this formula is m the normal isotropic index of refraction, Vo is the Verdet constant of the Faraday 
effect (we added a frequency factor u to be consistent with the experimental definition and to give it 
the proper dimension) which is described by the anti symmetric hermitean tensor $^=^6^5^. The 
third term is uniaxial birefringent and associated with the Cotton-Mouton effect [4]. The dielectric 
tensor £ = n 2 is given by, 

/ V 2 B 2 \ 2mV ( V 2 \ 

e (B) = m 2 + -2— I + + 2mM - + M 2 B 2 BB =el + + e M BB (2) 

\ CO 2 J CO \ CO 2 J 
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In typical optical experiments [11] at room temperature the Verdet constant is 90°/mm/T. In a 
field of 15 T the relative magneto-optical perturbation of the dielectric constant VqB /u ~ 0.002 for 
one single scatterer is still small. It seems therefore rather safe to ignore orders higher than B 3 in 
the dielectric tensor. Both the Faraday (order B) and the Cotton-Mouton effect (order B 2 ) remain 
present and no physics is lost. Note, however, that in multiple (and even resonant single ) scattering 
the magneto-optical effects may augment since they are basically proportional to the length of the 
traversed path. 

Two properties of the dielectric tensor are of fundamental importance for the rest of this paper. 
First, when all material parameters e, Vq and M are real-valued, the dielectric tensor is hermitean: 
£ (B) = e (B)* . This is crucial for multiple scattering since it will finally guarantee long-range 
diffusion. The second property is caused by the anti symmetric term only: £ (B) = £ (— B)* ^ £ (B)* . 
This property destroys time-reversal symmetry and modifies reciprocity relations. As such it is quite 
similar to the Lorentz force in electronic systems where the symmetry between electrons and holes 
will be broken. 

A point scatterer with dielectric constant £ (B) located at position r in vacuum is characterized 
by the hermitean operator, 

e(B,r) = I + <,[e(B)-l] |r) <r| , (3) 

where v can be interpreted as a typical physical volume associated with the point scatterer. The 
abstract Dirac notation |r) (r| is used for a local interaction with the profile of a delta-distribution 
S(r). The t-matrix of this particle is most conveniently obtained by noting that the Helmholtz 
equation for the electric field at frequency u resembles a Schrodinger equation with potential 
V (t,lu) = [I — e (B, r)] to 2 and energy to 2 . As a result the t-operator - here also a second rank tensor 
- is given by the Born series, 

t (B,u) = V (u) + V (u) -Go (u, p) • V (u) + • • • . 

Here G (cu,p) is the free Helmholtz Green's function G (cu,p) = l/[cu 2 — p 2 A p ], with (A p )^ = 
$ij ~ PiPj/p 2 - For a point scatterer the Born series can be transformed into an ordinary geometric 
series [21]. Using the matrix identity (|v| = 1), 

1 vv (I - vv) R - Q (ie ■ v) 



PI + Q (ie • v) + Rw P + R P 2 -Q 

it follows that 



u i 



t pp /(B,u;) = (t -t 3 )I + t 1 $ + t 2 BB , (4) 
in which t is the normal Rayleigh t-matrix with phase-shift a (u) , 

+ ( \ -4^rcu 2 6ir 

t (uj) = — 2 — exp(ia) sm a . (5) 



Uq — to 2 — |«rcu 3 to 



At low frequencies the phase shift goes to zero proportional to to . Near the resonance to the phase 
shift changes rapidly from — ir/2 towards ir/2. Up to orders B 2 the other t-matrices are given by, 



to 

U (u, B) = -—fit, 

D7T 



U ' - 2 U \ 2,3 



t 2 (u,B) = —(w-C)K-[— fi% 

D7T \D7T 



h(u,B) =-— (to-(^) ftl (6) 

D7T VD7T/ 
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For future use we define the dimensionless variables r/i = £;/ (— Im t ) and = ti/t . In Eq. (6) we 
introduced three other dimensionless variables which will be used in the rest of this paper, 



_ QitsfB 6ire 2 F B 2 _ _ QitsmB 2 

l) 2 uj 3 v ' (e — l) 3 uj 3 v ' (e — l) 2 uj 3 v 



The absolute value of the magnetic field has been included here and from now one we denote by 
B a normalized field with |B| = 1. The dimensionless variables /i, ( and w will finally determine 
the role of magneto-optical effects in multiple scattering. Their magnitude is discussed separately 
in Appendix C. The variable ji can be argued to be the dominant one, and essentially equal to 
the variable qpl ~ V e ffB£ (£ is some mean free path) introduced by Lenke and Maret [11]. Since 
the Faraday effect in our case occurs inside the particles, the effective medium Verdet constant is 
proportional to the number density of the particles. For that reason qpl is independent on the 
number density (in the dilute regime). 

The t-matrix does not depend on incoming and outgoing wavenumbers due to its point like 
origin. In this case we encounter anisotropy only in polarization indices and not in wavenumbers. In 
a magnetic field the resonant scattering matrix mimics the classical Lorentz model for the atomic, 
diamagnetic origin of the Faraday effect [11]. 

The Optical Theorem is an identity for the t-matrix that follows from energy conservation and is 
a mathematical consequence of e = e* . For our point object it is readily checked that, 



t (B,u) - t* (B,cu) = t (B,cu) • J2 AG o P) • t* (B. 



^t(B,cu)-t*(B,cu) 



In this paper we will frequently denote the anti-hermitean second-rank tensor A — A* by AA. The 

summation sign J2 1S a short notation for / d 3 p/(27r) 3 . 
p 

III. THE AVERAGE GREEN'S FUNCTION 

In this section we consider an infinite medium randomly filled (number density n) with the 
magneto-optical scatterers discussed in the previous section. The average field is determined by the 
average (Dyson) Green's function which - by translation symmetry - takes the form in momentum 
space, 



e(B,r)cu 2 -p 2 A p / pp , to 2 -p 2 A p - S(cu,p,B 



^ (P -G(p,BH(p-p') . (9) 



In the independent scattering approximation the self-energy S (cu, p, B) is given by nt pp (cu, B) . We 
note that the tensors S (cu, p, B) and A p do not commute. Using Eqs. (9) and (4) a rather lengthy 
tensor analysis reveals that, 



G(p,B) 



a 2 - 
with 



'q + S 1 *+ R(B ' p) + S(B ' p) 

+ 1 + detG-i(B,p) 



(10) 
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R (B, p) = - (S 2 - aS 2 ) (a 2 + ap 2 - S 2 ) BB+a (S 2 - aS 2 ) (B • p) (pB + Bp) 

-a 2 (a - S 2 ) pp - aSi (a - S 2 ) [($ p) p - p ($ p)] + S 2 (a - S 2 ) ($ p) ($ p) , 
S(B,p) =-S 1 (s 2 -aS 2 )(B-p)[B($.p)-($.p)B] , 

det G" 1 (B, p) = (a - S 2 ) (a 2 + ap 2 - S 2 ) - (B • p) 2 (s 2 - aS 2 ) . (II) 



We abbreviated 



a = u 2 — p 2 — (S — S 



and defined the various components of the self energy as S; = nti where the ti have been defined in 
Eq. (6) 

The complex dispersion law is determined by the poles of the Green's function, that is 

detG" 1 (B,p) = , (12) 

usually called Fresnel's equation [4]. With the explicit expression for the determinant given above 
we find two solutions in the upper sheet 1 . If we write z 2 = u 2 — S + S 3 and cos 9 = Bp, these are 
explicitly, 



2 2 
z -p ± 



S 2 + z 2 S 2 ) sin 2 9 ± V(S 2 + ^ 2 S 2 ) 2 sin 4 9 + 4S 2 (z 2 - S 2 cos 2 9) (z 2 cos 2 9 - S 2 

2(z 2 - S 2 cos 2 ^) 



= ^S 2 sin 2 9 ± y ^S 2 sin 4 5 + S 2 cos 2 5 + (n 2 
= ±Z X {B-p) + 0(B 2 ,n 2 ) . (13) 

In the third equality all birefringent terms are ignored and the propagating modes obey a dispersion 
relation similar as found in Ref. [7]. However, for directions nearly perpendicular to the magnetic 
field or for larger fields, birefringence cannot be ignored and the modes achieve an "ordinary" and 
"extra-ordinary" character. We emphasize that Si is a complex number and thus modifies both 
wavenumber k and scattering mean free path I according to p± = k± + i/2£±. In former treatments 
[6] [7] [8] the difference between £ + and £_ caused by the Faraday effect (even without Cotton-Mouton 
birefringence) has not been taken into account. This will cause many different results compared to 
previous work. We show the dispersion law for the low-frequency case, Fig. 1, and the resonant 
situation, Fig. 2. In the low-frequency regime the difference in wavenumber between the two modes 
is very small, and the Faraday rotation is almost entirely contained in the mean free path. As evident 
from Fig. 3 the accumulated phase shift between two successive collisions can nevertheless be large, 
because the mean free path is large at low frequencies. 

The following exact identity can be derived directly from the Dyson equation (9), 

G (p+|) - G* (p-|) = G (p+|) • [AS+2p q — pq — qp] G* (p-|) . 

Using the more compact Liouville notation (Appendix A) the becomes, 



frequencies for which z 2 — S 2 cos 2 6 = need special attention. They have one propagating solution given 
by a ~ S 2 / ' z 2 and one non-propagating solution z 2 = S 2 (z), with p along B. 
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0( a+ «^|)-0.(-«p-f)={G( a+ «^ f ) G .( u _« H )}. 

• [AS + 2p q - pq - qp - 2cufi + £>(nfi)] . (14) 

in which {AB}^ is considered as the four-rank tensor AijBik- Identity (14) has also been generalized 
for a finite frequency difference between the two Green's tensors on the lefthand side. This becomes 
necessary to discuss dynamics in multiple scattering. 

IV. THE AVERAGE INTENSITY 

The propagation of the average intensity - for electromagnetic waves more generally the complete 
set of nine Stokes correlations (^EiEj^j - is characterized by the Ladder diagrams (Fig. 4). In the 
Boltzmann approximation the Ladder diagrams are formally equivalent to the solution of the Bethe- 
Salpeter transport equation. If the t-matrices do not depend on incoming and outgoing momenta (as 
is true for our case), the Ladder diagrams for an infinite medium generate a simple geometric series 
in the four-rank tensor, 

Q ± (u, q, B, 0) = £ G (u + ^, p+|, B) G* - ^, p-|, ±b) • n t (u> + ^, b) t* (u> ° 



(15) 



Explicitly, 

Qf jkl (u, q, B, 0) = £ G ia (u> + ^, p+|, B) (G*) tk (u> - ^, p-|, ±B) 

nt sj (u + ^, B) (t*) lt (u - ^, ±B) . (16) 



For later use, the bottom line of the Ladder diagram is allowed to have reversed direction of the 
magnetic field. The Ladder diagrams sum up to, 

£± q, b, n> = „ t (* + 1, b) r (* - §, ±B ) • ^^^^ ■ (ir) 

Usually, this can be worked out by decomposing Q (0,q) = Q + 8Q (0,q) into the various eigen- 
functions and eigenvalues of Q thereby treating ^Q(fi,q) as a Rayleigh-Schrodinger perturba- 
tion [7] [22]. We want to note here that, due to the magnetic field, Q is not a normal opera- 
tor and does therefore not necessarily allow an orthogonal set of eigenfunctions. For this reason 
Rayleigh-Schrodinger perturbation theory cannot simply be applied here. In fact, the tensor Q is 
not necessarily diagonalizable, and at most isomorphic to a 9 X 9 Gauss-Jordan matrix. If we let 
AG = AG(r = 0) = E P AG (p) then Eqs. (8) and (14) imply that, 

Q+ • |AG) = |AG) , 

(AS| • Q+ = (AS| . (18) 

Since in a magnetic field AS ^ AG, this means that Qq has - on the basis of energy conservation 
- an eigenvalue 1 with geometric multiplicity 1 but with algebraic multiplicity 2. In our case the 
magneto-optical effects only occur in the particles. As a result the "effective medium" described by 
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AG is only influenced indirectly and AG^= (cu/37rz) 8ik + O (n). It is convenient to normalize both 
eigenfunctions and write 

|1)= ^ |AS) 

1 [(1 + Im 773) |I) — Im r]i |<fr) — Im i] 2 |BB)] , (19) 



A 

where rji was defined in the first section, and A = 3 (1 + 21m 773) +2 (Im rji) 2 — 21m r/ 2 ; ||C|| = (ClC) 1 / 2 
denotes the norm of the matrix C (Appendix A). 

The eigenvalue 1 of Qq implies that C + (q) will be singular in the hydrodynamic limit q — > 0. 
This is the origin of long range diffusion. Due to the modification of reciprocity relations in a 
magnetic field, Qq has no longer an eigenvalue 1 and - as it turns out later - neither an algebraic 
multiplicity 2. The tensor Qq is intimately related to Coherent Backscattering and is in that sense 
also an "observable". It will be discussed in section IV. 

In what follows we will focus on the diffusive 2D eigenspace of Qq and ignore 7 other dimensions 
in polarization space in which the propagation is non-diffusive. Let N + be the 2x2 Gram matrix 
of these vectors, constructed from the scalar product defined in Appendix A. For simplicity we first 
set = 0, and consider the perturbation in the variable q only. To this end we let 

+ 2 + (q) |1> (r\ Q+(q) \\)\ _^ J+ fq-D n -q «+ (B) + q • D 12 q 



P+(q)= N l * v ^ " * v ^ =N 

[q> \(\\ Q+(q) |r> <r| Q+ (q) |r> J ~ WD 21 -q q D 22 q ; ' 

in which 

k + (B) = (r| 1 - Q+ |1) ~ B 2 . (20) 
The Ladder tensor (17) can now be written as 

1 { (r\\ 



£ + (q) = ntt*- Mr), |1) 



N+-P + (q) \ (1| 



The matrix N + — P + (q = 0) is obviously in Gauss-Jordan form. After inverting the matrix 
N + — P + (q), the diffusive mode can be found by looking for the divergent entry as q — > . To 
associate a diffusion tensor to the diffusive mode, dynamics has to introduced in the form of a 
Laplace frequency as indicated in Fig. 4. (For = the diffusion tensor can never be obtained). 
The dynamical expression can be obtained without explicit calculation using the tensor identity (14) 
for O^O. The result is, 



£+ (q 0) - „tt* • |r> l - (II - ^l^L IMI 

(q ' )_ ' } -6^Q(||AS||y3y 1 + q-D 21 .q ( - 2 6 -.Q + q.D+(B).q 



(21) 

where the Boltzmann diffusion tensor is defined as, 

q • D+ (B) -q = 11 A ^ ^ q D 21 q = ±- (AS| SQ + (q, B) |I) . (22) 
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In the dynamic treatment we did not incorporate O-expansions of scattering matrices t(u ± 0/2). 
This becomes necessary if we want to be consistent with the correct conserved quantity [23] which 
is not (^EiE*^i as would follow from Eq. (21) but rather (eikEj-E*) . This will modify the transport 
velocity in the diffusion tensor, which is necessarily a scalar. Therefore the diffusion tensor is here 
basically (one-third times) the "transport mean free path tensor". Keeping in mind this aspect, we 
will nevertheless continue to use the more common word "diffusion tensor". Expression (22) can be 
shown to be equivalent to a DC Kubo-type formula for the "electromagnetic" conductivity tensor 
(i.e. density of states times diffusion tensor), applied to the special case of point scatterers. 

The lefthand eigenfunction | AS) must be a linear combination of the complete set of orthogonal 
eigenfunctions {|i)} - =1 9 of the hermitean tensor Q (B = 0, q = 0,0 = 0) obtained in Refs. [7] and 
[22]. A contribution of any of the last six of these will cause the diffusion tensor to be anisotropic. 
In fact is equivalent to |7) and |BB) is a linear combination of the first four (see also section 
IV.C.3.2). 

It is convenient to consider the symmetry of the tensor implied by reciprocity. Reciprocity of 
both wave and conjugate wave in the intensity requires the Ladder tensor C + to obey, in general, 

£± fc; (q,B,0) = £± fc (-q,-B,0) . (23) 

By applying this symmetry to the diffusion approximation found in (21) we conclude that the diffusive 
mode must satisfy 1(B) = 1(— B) and the diffusion tensor (B) = (— B) . The representation 
of the diffusive mode of £ (identified as the lefthand eigenfunction) must be a linear combination of 
the tensors I, <fr and BB , with real-valued coefficients. This will seen to be different for the diffusive 
mode of C~ . The general form of the diffusion tensor is given by, 

D+ (B) = Aso (B) I + (B) BB . (24) 

The coefficients are only dependent on the absolute value of B. We emphasize that the Ladder- 
diagram (21) only provides the symmetric part of the diffusion tensor. The anti- symmetric part 
is proportional to the tensor — i<f>ij=€ijkBk and responsible for a transverse Hall-type current, can 
be found only by looking directly at the current [16]. More generally, reciprocity guarantees the 

. r . i transpose 

"Onsager relation" D£ (B) = D£ (-B) which requires the anti-symmetric part to be odd 
in the magnetic field. 

A reciprocity relation between and can be obtained by transposing only the conjugate 
wave. From e (B) = e* (— B) it follows that, 

C^(p + p',B,0) = £^(q,B,0) , (25) 

where the tensor is given by the sum of the most-crossed diagrams (Fig. 4). Thus this reciprocity 
relation relates two sets of topologically different diagrams. 

The regime for which the diffusion approximation discussed here becomes valid can be character- 
ized by the criterion, 



det [N+-P+ (q,B 



K + B\ 



This is equivalent to q 2 <C 
of all other modes [7]. 



1 + 2 B-q 



£ 2 /5 . A similar criterion holds to guarantee the irrelevance 
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In the next subsections we will calculate explicitly the four-rank tensor Q^ 1 (cu, q, B) for our 
model and from that in particular the diffusion tensor. Next we translate the Ladder tensor to 
a slab geometry by adding boundary conditions. We will discuss polarization properties (Stokes 
parameters) due to the impact of the magnetic field on both the diffusive and the non-diffusive 
modes in transmission. 



A. Calculation of Ladder Tensor 

To find the explicit form of the four-rank tensor Q^ 1 (cu, q, B) defined in Eq. (15) one must perform 
momentum integrals over the product of two Green's tensors (9). The exact expression (10) is 
quite cumbersome for general multiple scattering calculations. In the Boltzmann approximation one 
considers only the "most divergent terms", that is the terms that become divergent for very low 
density. Keeping in mind that in all momentum integrals in the next sections the various divergences 
occur near the dispersion law (a 0) we can consider a ~ n . Furthermore S,- ~ n. Another 
dramatic simplification can be obtained by consequently expanding in B and B 2 . Altogether this 
replaces the Green's function by the expression, 

g (b ' p) iv + 5 - 1) Ap - § l (B 'P } - (§ - 1) h (B 'P } - 5 v (B 'P } ' (26) 

in which a = u 2 — p 2 — S and 

L(B, P ) = -$+($ p)p- P ($ p) , 
H(B,p) =BB+(B p) 2 pp-(B p)(Bp + pB) , 
V(B,p) =-(*.p)(*-p) , 

are three transverse tensors of rank two. Only the first is linear in the field and anti-symmetric, 
signifying the Faraday rotation of a spherical wave; H and V are birefringent parts quadratic in the 
field that preserve reciprocity. In the simplification procedure outlined above, the longitudinal part 
of the Green's function has been lost. As a result an ultraviolet singularity has disappeared and 
G (p — y oo) ~ 1/p 2 . Leading orders in the number density of the calculations in the next sections 
will not be affected. In what follows powers in the magnetic field higher than two will frequently be 
ignored without explicitly saying so. 

The radial momentum integrals for q = all take a form frequently encountered in speckle cal- 
culations [24], 

/ (n, m) = V / xm = (^ + ^- 2 ) ! » (-1)™" 1 (-t)"^- 1 . (27) 

V ' ; Y a o(«o) ir(n-iy.(m-l)l2 n + m u n + m - 2y J y J y J 

The q-expansion can be done using the radial momentum integral, 

V a o(P + q/2) Mp-q/2)] 

lPqJ tt (n-l)!(m-l)!2»+™+icu»+™-2 1 L) [ l) ■ [Zii) 

Here I = — u/Im S denotes the scattering mean free path in the absence of a magnetic field. Angular 
averages have to be performed for rank-two, rank-four and rank-six tensors. We find, 
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(PzPjPkPl) 
(PiPjPkPlPmPn) 



JE (fiijfikl + &ik&jl + &il&jk) = Vijkl 

7 i^ij^klmn ~\~ ^ik^jlmn ~\~ ^il^jkmn ~\~ ^im^jknl ~\~ ^in^jklm] 



(29) 



ijklmn 



Using this information the integrals can be carried out straightforwardly and we will only quote the 
final result, 

4tt 



A% kl (u>, B, q = 0) = — £ Ga (a;, p, B) G* lk (a;, p, ±B) 



Here, 



- nl - Im ^3) Tijki + j^m^jki t j^mFikji + \ 



2^ 



ii]2 ) ^jfc; 

(30) 



Tijkl 
J~ ijkl 

QijH 

Jijkl 

Sijkl 



(6SijS k i + 8 lk 8 3 i + SuSjk) , 
-QijSu - \ [$uS jk - 5 a $ jk - 5i k $ji + $,■*£,•/) , 

Y^^ij^kl 7 (Vijkl ~\~ B^B m ')^ rn j k l ~\~ B j B rn '\^ m ^ k l ~\~ B k B m ')^ rn ^jl ~\~ Bf B m ')^ rn ^j k ^ 

+ YE (^BiBjSki + BjBiS tk + B t B k Sji + B t BiS jk + BjB k S t i) , 
± (SijSki - BiBjSki) + ^ ($jfe,-$/j + $/,-$jfej) , 

+Tg (SjkBiBi + S t iBjB k — 8jiB t B k — 8 lk B 3 Bi) , 



It is convenient to make an isomorphic transformation to the 9-dimensional vector space of 3 X 3 
matrices. We will restrict ourselves to the three-dimensional subspace spanned by the hermitean 
matrices {I, BB}, in which the long range diffusion occurs. The tensor Q ± (uj,c{ 
closed in this subspace. This base is not orthonormal and has Gram matrix, 



0,B) 



is 



M 



/3 1 \ 

2 

1 1 



With respect to this base the tensor A^ for q = takes the form, 



4±(cu,B,q = 0) 



i|.4(u;,B,q = 0) \j) 



]TG(cu,p,B)G* (u,p,±B)|j) 



(31) 



with, 



Af 1 = 2- |Re nl - 21m r, 3 + flm r] 2 ±\ \ Vl f 



A ± 

1± 

l 22 

1± 

l 23 

lf 3 

1± 

1 33 



A d 
A d 



2 _ 

3 5 
l± 



I Re 77^ + ^Im ?7 2 - flm rj 3 ± | | 

32 - "T^fal T»7l) 

■31 - i - i5" Re ^i 2 + il Im ^ - |im ^3 ± ^ M 2 

8 ±,Re V l + ilm ^ 2 8 



A : 

A. 



15 



35 



15 



3 

Im r/3 



This matrix is symmetric and for the "+" choice even hermitean. We will also need the four-rank 
tensor nt (B) t (±B)* with respect to this base. If we denote, 



12 



6ir 



Uh(u,B) = — (i\t(u,B)t* 



(",±B)|j> , 

we obtain the symmetric (and for the "+" choice again hermitean) matrix, 



T± 
33 



3 - 6Re 6 + 2Re 6 ± 2 |6| 2 Uf 2 = U 2 \ =2 fc± 



2-4Re ^ 3 ±2|^r 
1 - 2Re £, + 2Re f 2 



'12 

T± 
13 



u£ 2 

T± 
31 



With respect to this base the Ladder sum becomes, 



£± (w, B, q) = nt (B) t (±B)* • Ml) , |$) , |BB) 





1 - 2Re £, + 2Re f 2 



( (II \ 

(bbi; 



l± (q) 



where (q) 



M 



A± (q) • M" 1 • . For q = we find, 



L 



li 



r± 

-^12 

r± 

-^21 







22 



-^23 

r± 

-^32 

r± 



13 



31 



L 



33 



| Re + 31m r] 3 - Im ^ =F f | 2 
+ 6Re 6 T 2 |6| 2 " 2Re 6 + f (liT^) (f i ± G 
|* ('7iT^7i)-2(ei±? 1 ) 

H^t^)- (6±ei) ±o 

1 + ^Re ^ 2 - flm r] 2 + Im ^ 3 T ^ ki| 2 

+ 2R e e 3 Tieii 2 + !^iT^)(ei±ei) 

1 '(»7i T»h) 

(^T^)-|(ei±ei) 

Re ri\ - flm ^ 2 + Im^T^ | 2 - 2Re 6 + 2Re 6 
10 Re T]l - | Im ?? 2 + Im r] 3 =p ^ ^ | 2 

+2Re 6 T | 16 1 2 " f Re 6 + ^ (»7i T ^i) ± li) = 
| + |jRe ^ - §f Im r] 2 + flm ^3 - f Re 6 + f Re 6 



10* 

X 

10 
X 
10 



(32) 



(33) 



This matrix is not symmetric. Some entries as well as the determinant for the "+" choice are zero 
on the basis of the Optical Theorem (8). 

An extensive calculation - of which we shall not mention any intermediate results and simply 
quote the final result - yields the q-dependent part of (q) with respect to the same base 



£L± (q) = f 2 [S V + T ± (B • q) 



(34) 



in which the isotropic part reads 
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St = 1 - 3Re r,l - 31m f] 3 + flm r, 2 ± | |^ | 2 

-2Re 6 ± | 16 1 2 + | Re 6 - ^ (liT%) (6 ± 6) 

^2 = -f- i(viTVi) + lki±li) 

Sti = -^(^iT^) + |(ei±ei) 

51 = | - f Re vf + flm r ]2 -llmr ]3 ± 1 -i \ V] \ 2 

-§Re6 ±|ieii 2 -^(^T^)(ei±ei) 

St 3 = | - fRe vl ~ flm + flm V2 ± £ l^ 2 - |Re 6 + fRe 6 
5 3 1 = | - fRe - flm ^3 + flm m ± £ l^) 2 

-fRe 6 ± i 16 1 2 + §Re 6 - T^g (»7i T li) (6 ± £i) 
^3 = i " iRe + flm V2 - flm Vs + §Re 6 - §Re 6 

and the anisotropic part, 

T± = -|lm ??2 ±|| ??1 | 2 ±||6| 2 -fRe6-|^(^iT^)(ei±ei) 

t& = -^(^T^) + |(ei±ei) 
r 2 i = -fi^iT^ + f (ei±ei) 

T± = f - |Re r? 2 + ilm V2 - flm ^3 ± f k_i| 2 
-fRe S ± f 16 1 2 - f§^ (^7i T ^x) (6 ± ^1) 

t 3 % = -^(^T^) + i(6±ei) 

Ti| = -f - iRe rft + flm r, 3 - fflm r, 2 ± i M 2 + fRe £3 - fRe £2 
Tal = -f - iRe ?? 2 + flm ^ - fflm ^ 2 ± i M 2 

+f Re 6 ± i 16 1 2 - iRe 6 " & foi T Vi) (6 ± 6) 
T & = "i + § Im ^ - §Im ^2 - iRe 6 + iRe 6 

Finally we give the perturbation of the tensor Q^ 1 in the dynamic variable 0. Using G (to ± 0/2) = 
G (cu) ± cufi dG/dcu 2 , the most divergent terms (proportional to 1/n) can be obtained straightfor- 
wardly using the standard integral /(n, m) introduced earlier in this section. If we represent the 
matrix elements of SQ^ 1 (0, B) with respect to our base by, 

SL* (0) = -zO£W ± (B) (35) 

we obtain, 
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3 — | Re r\\ — 61m r/ 3 + 2 Im r\ 2 ± | |??i| 2 

- 6Re 6 ± 2 |6 1 2 + 2Re 6 -i(^T^) (6 ± 6) 
-i(»/iT»/i) + 2(ei±? 1 ) 

-»(liT^)+ (&±£i) 
1 - ^Re ri\ + flm r] 2 - 2 Im r] 3 ± ^ ^ | 2 

- 2Re 6 ± I6I 2 - ^(/7i T^7i) (6 ±ei) 

1 - ^Re r? 2 + flm r] 2 - 21m r? 3 ± ^ ki | 2 + 2Re 6 - 2Re 6 
1 - ^Re T]\ + flm r] 2 - 21m ?? 3 ± ^ | 2 
-2Re 6 ± | 16 1 2 + fRe 6 - f (liT^) (6 ± £i) 
| - |jRe ^ + flm ^ - flm ^3 + fRe 6 - fRe 6 

The identity (AS| 8Q + ( fi,B) |I) = — 6cuO + 0{n) follows from the Optical Theorem and was 
used earlier for the calculation of the Boltzmann diffusion constant. The explicit calculations (35) 
have been verified to obey this identity. Since we ignored terms higher order in density in this 
calculation, the influence of the scatterers on the transport velocity of multiply scattered light [23] 
in a magnetic field has been neglected. 

B. Boltzmann Diffusion Tensor 

With the calculations in the previous subsection and using the eigenfunctions in Eq. (19) we find 
that the Boltzmann diffusion tensor (22) is, 

q • D+ (B) -q = ^ [(1 + Im r, 3 ) SL^ (q) - Im ^1+ (q) - Im ^1+ (q)] , 

with the explicit result that the symmetric part is given by, 

D+' sym (B) = U [d iso I + 4„iBB ] , (36) 

d iso = 1 - 3Re r\\ - 21m r] 3 + flm r] 2 + f \r] x | 2 - f (Im r] x ) 2 
-2Re & + f |6 1 2 + fRe 6 + flm ^Re 6 , 
4ni = -flm V2 + f l^l 2 - f (Im ^) 2 + f |6 1 2 - fRe 6 + flm ??iRe 6 • 

In Appendix C it is shown that the variable ji is the most important one. Ignoring ( and w in 
Eq. (36) yields that d- lso — 1 ~ c/ an ; ~ /i 2 ~ B 2 . The dephasing in the most-crossed diagrams will 
later be shown to be determined by the same variable ji. In our model both phenomena are thus 
seen to be of the same order of magnitude. This dephasing has been measured already [10]. For the 
experimental verification of anisotropy in the diffusion tensor of light in Faraday-active media we 
may thus be optimistic. 

In Fig. 5 we show the modifications to the symmetric diffusion tensor as a function of the phase- 
shift a of the Rayleigh t-matrix defined in Eq. (5). At low frequencies and at resonance it is inferred 
that the diffusion constant is suppressed, and relatively more perpendicular to the field direction. 
At low frequencies (a ~ 0) - a regime where our model should coincide with the exact Mie solution 



W? 2 

w 2 \ 



22 



W ± 



23 



w. 



w. 



32 

± 
13 



31 



W ± 



33 
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at low frequencies - we find AD|j~/A.D^ = 1/2. At resonance (a ~ 7r/2) this is changed into 1/3. 
For the Beenakker-Senftleben effect one also finds the modification to be largest perpendicular to 
the field, and for the heat conductivity - at low fields - a ratio AA||/AAj_ = 2/9 [15] is obtained (for 
large fields this saturizes to 2/3). At intermediate frequencies we infer that the diffusion may in fact 
increase in a magnetic field. For the Beenakker-Senftleben effect one always finds a decrease. 

C. Light Propagation in a Slab Geometry 

In this subsection we calculate multiple light scattering in an optically thick plan-parallel slab im- 
posed to a homogeneous external magnetic field. This is the simplest theoretical model that includes 
all physics and is experimentally relevant. The angular transmission tensor of a slab with thickness 
L in (nearly normal [25]) direction k and normal incidence is given by [26] (y = (L — z 2 ) k, z = zizj 

T z , k = f L dz x f L dz 2 (Airy) 2 G (y) G* (y) ■C(L-z 2 ,z 1 ) ■ (Attz,) 2 G (z) G* (z) . (37) 

JO JO 

Similarly the reflection tensor is, (y = z 2 k, z = zizj 

n z , k = f L dz x f L dz 2 (Airy) 2 G (y) G* (y) • C (z 2 , z x ) ■ (Attz,) 2 G (z) G* (z) . (38) 
Jo Jo 

Here C(z 2} Zi) is the Ladder tensor for the slab geometry (integrated over transverse coordinates) 
and must be found by adding boundary conditions to the bulk results obtained earlier. The Green's 
tensors G in the equations above are Dyson Green's tensors and the Fourier transform of G(p) defined 
in Eq. (9). They decay exponentially in space and thus signify the light propagation in a narrow 
"skin" layer, that is the first scattering mean free path for the incident wave, and the last one for the 
emerging wave. For the Ladder tensor in the bulk we obtained in Eq. (21) £ (q) ~ |d) p (q) (d| with 
p (q) essentially the well-documented scalar Ladder sum but with an anisotropic diffusion tensor. 

The boundary effects caused by the finiteness of the scattering geometry are coded in the 
Schwarzschild-Milne integral equation for C(z\, z 2 ) that can be derived from the Bethe-Salpeter trans- 
port equation. For scalar waves and isotropic scattering this was investigated by Nieuwenhuizen and 
Luck [17]. The generalization to anisotropic scattering was recently carried out by Amic, Luck and 
Nieuwenhuizen [28]. These methods are not yet applicable to vector waves in anisotropic media. 
Therefore we will use the (improved) diffusion approximation in which the boundary conditions are 
introduced by the imaging method [7] [26]. This consists of adding trapping planes for p(z } z') at 
some distance z beyond the slab boundaries, called the extrapolation length. For isotropic media 
with k-anisotropic scattering it known from Milne theory [29] that z 0.710£*, that means always 
essentially one transport mean free path. A similar conclusion is reached by considering the bulk 
diffusion equation subject to the "radiative boundary conditions" that the incident incoherent flux 
on both sides of the slab boundary vanishes (see e.g. Refs. [30] and [31]). This method can also be 
formulated for vector waves with an anisotropic diffusion tensor in a slab geometry, 

(— V • D • V + a 2 ) p (z, z') = source, 
\p(z,z')-±z-I)-Vp(z,z') = z = } (39) 
\p {z,z') + \z-T) ■ Vp(z,z') = z = L, 

The source is usually supplied by single scattering; a 2 is an absorption rate to be specified later. In 
this approach the boundary conditions effectively locate trapping planes at distances z = 2D ZZ = 
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-^£* zz outside the slab boundaries. The relevant diffusion constant for the slab geometry is thus D z 



[27]. 

The solution of Eq. (39) can be found in literature. For our purposes we need p (z, z' 
which reads, [26] 



for 



a 



p(z,z') 



1 (L + 2z 



j\\ 2 



j\ 2 



4D, 



L + 2z 



(40) 



The angular incoherent reflection coefficient for L = oo (at approximately backscattering) adopting 
first ordinary scalar skin layers can be found straightforwardly by integration of Eq. (38) with scalar 
Green's functions. The same holds true for the incoherent transmission (at forward scattering) [32] 
[33], which we shall need later also for a ^ 0. We have, 



R (back) 
T (forw, L) 



£ 1 + x — exp (—2xz /£) 
2D~7 Z x(1 + x) 2 

1 (z + £) 2 L/L a 
D zz L + 2z sinh L/L a 



a£ 



L 



oo 



a 



'D, 



(41) 



where z ~ 2.1 D zz . 

We still must deal with the polarization anisotropy caused by first and last skin layer. For that 
we need the Dyson Green's tensor in real space G (r) . Since the skin layer is roughly one mean free 
path in thickness and k£ ^> 1, the far field suffices. If we restrict ourselves again to low fields and 
low density we need the Fourier Transform of expression (26), in which p will be transformed to r 
and the remaining integrals over the absolute value of p are standard. We find, 



Gfr,Bl 



exp (ikr — r/2£) 



1 . v 
, 2 ' £ 



Airr 

r l 9 /r\ 2 \ A 1. r T __ ( 1 
V M- A P + -ini-L(r.B) + - 

I 8 11 \£J J 2 1 £ y ' ; V 8 



1 . r 
7T^2 



£ 



H(r,B) 



V(r,B) 



(42) 



The four transverse tensors herein have been defined in section IV. A. 

One result can be obtained without going into deep calculations and concerns the "all channel 
in all channel out" incoherent transmission coefficient T . This is the total transmitted energy flux 
integrated over all incident and outgoing directions. This transport quantity is intimately related to 
the conductance for electrons by means of the Landauer formula [34]. It can best derived from the 
diffusion formulas (39). An even simpler picture is obtained by replacing the source by a boundary 
condition as well. If u (z) is the radiation density then, 



-d 2 z u (z) 

Jz (m,0) 
Jz (in, L) 
J z (out,L) 

^From this it follows easily that, 





b 

T(L 



(o)-i^V(o) 



L) + \Dt z u< 



[L) = 



(43) 



T fB, L) 



4Dfz (B) 
L + 4D+ (B) 



(44) 
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This formula also comes out in exact Milne calculations for angular anisotropic scattering without 
the diffusion approximation [29] [17], though with a somewhat different numerical value for the 
extrapolation length. On the basis of this formula the all-channel transmission depends on both 
direction and absolute value of the magnetic field. This is a clear experimentally verifiable prediction. 
In fact, the all-channel transmission is the only known transport quantity that depends solely on (the 
anisotropy of) the diffusion constant. In the presence of absorption the diffusive absorption length 
L a ~ Df z is also dependent on magnitude and direction of B. In Fig. 6 we have displayed 
the difference in transmission for two perpendicular directions of the magnetic field. The difference 
Tj| — Tj_ can have both signs, depending of the phase shift. At low frequencies (zero phase shift) the 
total transmittance decreases as a function of the magnetic field, for all directions of the magnetic 
field. This is similar to normal electronic conductors. Note however in Fig. 5 that there exists a 
phase-shift region where the diffusion constant and thus the longitudinal conductance increase in a 
magnetic field. 

The complication in reflection is that, contrary to transmission, short paths also contribute so 
that in principle the full Ladder tensor is required, and not only its diffusive approximation. As a 
result polarization effects occur even without a magnetic field. A realistic treatment of polarization 
effects in reflection requires to go beyond the diffusion approximation. In the subsections that follow 
we therefore only discuss the transmission. 



1. Stokes parameters in Diffuse Transmission 



For light incident along the z-axis, we can use Eq. (42) to transform the four-rank tensors GG* 
in Eq. (37) into a transfer tensor M^ of the skin layers acting as, 

T z , k = Mt • dz x £ dz 2 e~ T -l C{L- z 2 , z x ) e" Tl j • M+ . (45) 

With the boundary conditions on the diffusive bulk outcome (21), the integral gives the familiar 
scalar result mentioned in Eq. (41) so that, in forward direction, 

T k , k ~ \ x Mi • |d> <d| • Mi . (46) 

where the diffusive eigenfunction |d) = \/A |1) was found in Eq. (19), and Zq ~ 2.1Df z . The skin 
layer polarization tensor in direction k is given by, 



M 



± 



^f 2 Re r\\ - /ilm n 3 



A k A k + ^i/i [?7iL k A k =p ?h A kL k ] + [^A k H k - ?? 2 H k A k ] 
+ [r7 2 A k H k - r? 2 H k A k ] - -J 2 [^Vk A k + r\\ A k V k ] ± \f 2 |^| 2 L k L k . (47) 



We allowed - for future need in the most-crossed diagrams - a possible opposite sign for the magnetic 
field in the complex conjugate wave. The scalar constants /,• arise from the space-integrals of the 
various factors (r/£)' in Eq. (42), and depend on the direction k. For details we refer to Appendix 
B. We want to emphasize here that by arbitrarily setting /,• = there will be no skin layer effects on 
polarization and the result Mi = A k A k is obtained, as in Ref. [7]. On the other hand putting /,• = 1 
the skin layer behaves exactly as the bulk in the sense that (r) = I and (r 2 ) = 2£ 2 [36]. Using the 
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identity (14) this would imply that the skin layers transform both left eigenfunctions |d) in Eq. (46) 
into the transverse tensor J2 p AG (pk)^> ~ l^k) + O (n) . In that case we would have to conclude 
that the polarization obtained in the bulk due to magneto-optical effects will be compensated exactly 
in the skin layers, and that we finally end up with no polarization as without a magnetic field. From 
this discussion it may be evident that the anomalous step length distribution in the skin layers plays 
a crucial role for the emerging polarization. The variables /,• are fixed by the boundary conditions. 
In Appendix B we calculate that fx = 1.58 and f 2 = 2.16 in forward transmission. This brings us to 
the interesting conclusion that polarization will be present in diffuse transmission. 

Polarization is best described in terms of the Stokes parameters /, Q } U and V which we can 
obtain from Eq. (46). Since Eq. (46) separates the incoming and outgoing Stokes variables completely, 
working out |EE) trans ~ Ml^ • |d) will suffice. We obtain, 

T k , k ~ Mt ■ |d) ~ (1 + Tx) | A k ) + T 2 |L k ) + T 3 |H k ) + T 4 |V k ) , (48) 

where 

Tx = -Im V3 (fx - 1) + (Im ^ f (f 2 - fx) ; T 2 = -Im ^ (fx - 1) 

T 3 = +Im f] 2 (fx - 1) - (Im m ) 2 (/ 2 " fi) ; T 4 = (Im m ) 2 (f 2 -fx) 1 } 

The Stokes parameters can now be found by choosing a frame with z-axis along k and x-axis and 
y-axis perpendicular. Then T xx = Tx + T 3 B 2 — T4B 2 , T yy = Tx + T 3 B 2 — T4B 2 and T xy = 
(T 3 + T 4 ) B x B y - iT 2 B z , T yx = (T 3 + T 4 ) B x B y + iT 2 B z . Hence, 

= T 2 B Z = -Im r,x (fx - 1) (B-k) , 

(T 3 + T 4 ) (B 2 X - B 2 y ) = Im V2 (fx - 1) (l - (B-k) 2 ) cos2</> , 



T _ ~T + T 

1 - 1 - xx \ - 1 - yy 

T -T 

xx - 1 - yy 

1 ~ T + T ~ 



U T xy + T yx ( ^ + ^ = Im ^ ^ (B-k)') sin 24> . (50) 



/ T +T 



yy 



The Stokes parameters V, Q and U do not depend on the particles density, and only on the direction 
of observation with respect to the magnetic field and the frequency. A circular Stokes parameter 
V ~ B persists when the magnetic field is along the slab. A linear polarization proportional to B 2 
remains if the magnetic field is perpendicular to the slab. The angle (j> locates the direction of the 
magnetic field in the xy-plane. Without a magnetic field the Stokes variables would decay to zero 
exponentially as exp ( — L/£**) [8], where the depolarization length £** is associated with non-diffusive 
eigenvalues of the Ladder tensor. We emphasize that the imaginary part of the variables r/i signifies 
the tensor nature of scattering mean free path and scattering cross-section . If this effect were not 
taken into account, the Stokes parameters would decay exponentially to zero in a magnetic field 
as well. In leading (zero) order of the particle density, the polarization in diffuse transmission is a 
surface effect. 

Finally we want to note that the Poynting vector associated with the transmission tensor (48) 
can easily be shown to be equal to J k = k Tr T kjk . If we assume that the incident light is completely 
unpolarized, an additional similar trace factor comes in due to the incident skin layer. In forward 
transmission we thus obtain for the Stokes parameter /, 



T (forward) 



4 + 

D zz 



1 + Tx + l - (T 3 + T 4 ) (l - (B-z) 2 



2 



i • < 51 » 
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This differs from the all-channel transmission (44) but nevertheless depends on the direction of the 
magnetic field. In Fig. 6 we show (T\\ — T±^j /T for both the all-channel transmission (44) and the 
forward transmission (51), quadratic in the field. 



2. Stokes Parameters in Non-diffuse Transmission 



In this subsection we discuss the impact of a magnetic field on polarization modes that are 
not subject to long-range diffusion. These modes are known to carry polarization characterized by 
either a non-vanishing U, Q or V Stokes parameter. In the absence of a magnetic field they decay 
exponentially in space, the characteristic length scale being associated with non-diffuse eigenvalues 
of the collision operator [22]. Martinez and Maynard [8] calculated the Miiller matrix for Rayleigh 
scatterers after a given number of scattering events n using a transfer matrix method. They showed 
that the magnetic field pronounces the exponential decay of polarization in order B 2 and in addition 
causes a rotation proportional to VoBn in the (Q } U) plane. Nothing was seen to happen to the 
circular polarization V. 

We will now investigate the decay of polarization in our transport theory. Contrary to the diffusive 
mode, terms linear in the field prevail. Therefore we restrict to modifications linear in the field. In 
this approximation the tensor Q^ 1 (q = 0, B) reads, 

Qo ( B ) = 7^ (MijSki + + SuSjk) + j^iriiFijki T j^ifiiFikji + ^i^S ± Tfi^ Z iSi • ( 52 ) 
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Here Ti 3 u has been introduced earlier in Eq. (30). Furthermore 

z\f M = 6$ij5 k i + §ij5 lk + 8 a $ kj ; z\f kl = 6$^^,- + $ij8 ik + 8 kj $u . 

For B = the tensor 1 — Qq (B) is hermitean. Its nine eigenvalues and orthonormal eigenfunctions 
associated with the various Stokes variables have been found in earlier work [22], 



Ai = |1>* = W3 (/) 

A 2 , 3 = 3/10 |2, 3) tk = 8 ik e^ k -^ 3 /V3 (Q) 

A 4 , 5 ,6 = 3/10 |4, 5, 6) tk = [8 ia 8 kh + 8 ib 8 ka ] /y/2 (U) 

X 7 ,s,9 = 1/2 |7, 8, 9) tk = i [8 ia 8 kh - 8 lh 8 ka ] /y/2 (V) 



(53) 



For the last six eigenvectors we have (a, b) = (1,2), (1,3) and (2,3) respectively. We shall take the 
direction of the magnetic field along the 3-axis. With respect to the base (53) the tensor 1 — Qq (B), 
occurring in expression (17) of the Ladder diagrams, takes the form, 
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in which \i = (3 + /2\/6 and X2 = (l + i\/?>) /2y/6. We introduced 



7 (6 ± £ x 
5 6 ± £i 



z ° = 7e . z =F - 2 (fi ± f J_ ^3 = i {m T rji, 
Zl= 7e[ l =F »7i) — 4 ± f J] ^ 4 = i {m T rji] 
z 2 = 2{r ]l ±r] 1 ) + 7i (6 =F ^5 = i(6T ^ 

For the normal Ladder diagrams all z,- are real- valued and z = 0. We shall focus on two 
parts of this matrix outside the one-dimensional diffusive domain (given by some linear combi- 
nation of first and seventh eigenvector). Firstly the "linear polarization subspace" spanned by the 
{|2) , |3) , |4) , |5) , |6)} . Secondly the circular polarization space spanned by |8) and |9) . For brevity, 
we do not discuss the (birefringent) coupling between both, transforming linear and circular polar- 
ization in general into elliptically polarized light. 

The matrix L (B) restricted to the five dimensional non-diffusive subspace with linear polariza- 
tion can be diagonalized. It is inferred that the magnetic field lifts the five-fold degeneracy of the 
eigenvalue 3/10. We find for the eigenvalues and eigenvectors, 

/ 1 \ 



A 2 = — 
10 



e 2 
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72 
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+ §n/3 
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3 4- 1 ' 
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/ ±iXi \ 
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( ° ^ 



e 5 ,e = ^ 

\±ij 

Quite conveniently this set is again orthonormal, although the eigenvalues have become complex- 
valued. As a result, the change of the eigenvalues for a small variable q / can be obtained with 
Rayleigh-Schrodinger perturbation theory, 

a,-(q,B) = -(\ i \5Q ± (< l ,B)\\ i ) . (54) 

For B = most matrix elements in Eq. (54) have already been obtained in Ref. [7]. Up to terms 
linear in the field and for low density we have, 



-<^(q,B) 



p q) 



A tJ A k l + ^Aim^mjAkl ± ^Aij^lmAr 



-i (rjiLijAu =F T] 1 A lJ Li kj 



(55) 



As before A,-j 



Pifij and Lij = — <&i 3 + &i n p n Pj — QjnPiPn are two transverse tensors used earlier. 
The angular average can be expressed into the tensors Vijki and Wijkimn introduced in Eq. (29). The 
matrix elements (54) have been obtained with a symbolic processor. We find, 

a 2 (q,B)/£ 2 = ijq 2 -^ 

a 3 ( q ,B)/£ 2 = i§q 2 + ^3 2 + (eiTe 1 )(Mq 2 + HV^^±^)(ilq 2 + ^3 2 ) 

(56) 



a 4 (q,B)/£ 2 = ±§q 2 + j^ 6T(, gq 2 + a- 8 (^± i ? 1 )(fiq 2 + ^; 



a 5 (q,B)/£ 2 



70' 
35^ 



L 4 qi + Ui=f£i. 
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J-o 2 ^ 

28^3 J 



q T 35 y 3 , 



,56 



SXe (q, B) /f = A q 2 _ 1. 2 _ U T Q A q 2 _ _L ? 2 - l (r ]l ±rj 1 ) |q 2 + § g 2 
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For the Ladder diagrams we find that only A 2 (B) is real-valued and independent of B. The other 
four occur in complex conjugates. The eigenvalues of Q~ (q, B) remain real-valued. 
Writing A,- (q, B) = A,- (B) + Di (q, B) q 2 and by using the formal identity, 



dq 



exp (iq ■ r) 



■ exp 



L,(r,B) 



D t (q) q 2 + A, 4tt A (?) r r \ V D i (?) , 
the Ladder sum for linear polarization can be transformed into real space. The result is, 



(57) 



( U (r,Bl 



L P (r) ~ U- 







Ls (r,B) 



U (r,B) 



L 5 (r,Bl 



V 



(58) 



in which the unitary matrix U is given by, 



/ 1 i X i -i X i \ 
-\ + iXi -iXi 



U =~r 
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1 1 

i —i ) 



The complex value of the eigenvalues A,- (B) and "diffusion constants" Di (r, B) will give rise to 
oscillations in the Stokes parameters Q and U as a function of B and r. To determine these we can 
consider the special case that initial and final wave vectors are directed both along the B— axis and 
find Q (r) and U (r) . Although this is still a calculation valid in the bulk, it will approximately be 
valid for transmission in a slab geometry with length r = Lz. If the initial light is 100 % linearly 



polarized (U = 0) we can take the initial polarization matrix 
Using Eq. (58) we get for the propagated light, 



E i (°) E 3 i°] 



Ei (r) E 3 (r; 



2i 

^From this expression one can infer that, 

Q(r) 



+ {(^^)M r » + (^^)^±^},3),, 

|4>„ ■ 



+v ^£ j( r)- M r) 



I(r 
U(r 



~ e 



~ r/lp cos (k P r + tfo) 



~ e r/ ^sin(A;pr + <p ) . 
l(r) 



(59) 



The linear Stokes parameters are exponentially small in the optically thick regime as is well known 
from other work. In a magnetic field the U and Q acquire an oscillatory character: the linear 
polarization rotates along B proportional over the traversed distance. Although this may be evident 
in view of the Faraday effect standing at the base of this, we nevertheless want to emphasize that we 
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are addressing a (damped) diffusive mode and not the coherent beam. Since incident and outgoing 
wave vector are both parallel to B, the diffusive contribution (50) vanishes, and Eq. (59) contains 
the only linear polarization of the light. 

^From Eq. (56) we find the explicit result that lp and kp, determining the propagation of the 
linear polarization, are given by 

/23 1 /21 /- 133^ 10 T ,\ \i 1 21 /133 112 . 2 



£ p = h— ■ kp = -\—{ Rem Imfi = ^W— sin 2 a . (60) 

A/ 21 ' £ V 23 V 138 ' 3 U / £ V 23 V138 23 ) y ' 

In the low frequency regime the variable kp does not depend on frequency and has the value kp = 
1.29 X 9y/e/ (e + 2) 2 fVoB. This is 1.29 larger than the rotation of linear polarization of the coherent 
beam. 

A similar calculation is possible in the two-dimensional subspace spanned by the eigenvectors |8) 
and 1 9) which carry circular polarization. 

M,, b> = i + ^ + f & - ±4 - (ft * f ~ ^ " ■ <* ± ?.) (^ + } • 
a, (,. b) = I - i.,, + e {i q * - ±i + ({, t 1 ,) (i q 2 - ^) + > ± * ) + ^d) } ■ 

For the normal Ladder diagrams these eigenvalues are again each other complex-conjugate. Writing 
A 8 (q, B) = A 8 (B) + Dg (q, B) q 2 and L 8 (r) as in Eq. (57) we obtain in real space, 

The eigenvectors |8) and |9) couple in when the magnetic field is perpendicular to the direction of 
propagation. By taking a circularly wave as initial wave it is easily shown that the Stokes variable 
V (r) is given by 



Y^l ~ e ~ r/lc cos (k c r + ip ) , (64) 



Explicitly, 



e c =ej- ; k c = \-\ -Rem = -~\ - (l - 2 sin 2 a 
V5 ' £8V2 ' £ 8 V 2 v 



We thus infer that the amount of circular polarization oscillates as well as a function of distance. 
However, for the circular polarization the exponential decay is much faster and the wavenumber 
associated with the oscillation is much smaller. 



V. COHERENT BACKSCATTERING 

After the theoretical predictions in Refs. [6] and [7] the Coherent Backscattering in a magnetic field 
has been studied intensively by Erbacher, Lenke and Maret and a vast amount of experimental data 
is now becoming available [10] [11]. So far the cone is the only multiple-light-scattering phenomenon 
investigated thoroughly in relation with a magnetic field. As has been predicted theoretically the 
enhancement of the cone was seen to go down proportional to the magnetic field and good agreement 
was found between experimental data, theoretical predictions and numerical simulations [8]. 
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New features emerged from the experiments which mainly concern very high fields, up to 30 Tesla. 
The decrease of the Coherent Backscattering Peak due to the Faraday effect turned out to be very 
well described by the exponential step length distribution between successive collisions [36]. So far 
Coherent Backscattering in a magnetic field seems to be the only phenomenon explicitly sensitive 
to the whole distribution and not to only its two lowest moments [38]. In addition, a saturation of 
the enhancement factor was seen to occur for very large fields that also comes out of the numerical 
simulations by Martinez and Maynard. It turned out that only one single dimensionless variable is 
needed to describe the dephasing up to high fields, being the product of an average Verdet constant, 
the magnetic field strength and some mean free path. This is quite similar to the Beenakker- 
Senftleben effect for which the magnetic impact on the heat conductivity turned out to be a function 
of only the product i? = u p t of precession velocity of the magnetic moment and the molecule's mean 
free collision time [12] [14]. It seemed necessary to introduce a length scale £* F different from both 
scattering and transport mean free path. It describes the loss of circular polarization in multiple 
scattering (being relevant for the average Faraday rotation between collisions). Due to spin-flips it is 
in general larger than the transport mean free path (which determines the memory of momentum). 
Initially some optimism existed to have the Faraday effect enhanced by resonant scattering [10] but 
restoration of a numerical factor in the formulas made the effect less pronounced, leaving it as an 
open question. Our approach will contain a first theoretical attempt to incorporate the effect of 
resonant scattering on the Faraday rotation of polarization. 

Concerning polarization it is found that in opposite helicity hardly no Coherent Backscattering is 
found (even for large magnetic fields). This is in agreement with predictions first made by MacKintosh 
and John [7]. Furthermore, the azimuthal line shape of the cone was seen to remain perfectly 
rotationally symmetric. For linear polarization channels there is a special angle between incoming and 
outgoing polarization vectors for which the cone is maximal. This has been explained by the Faraday 
rotation in the skin layers. In all these experiments the magnetic field was directed perpendicular to 
the slab, which is the most convenient set-up in the magnet. 

To interpret the experimental data one focuses only on the diffusive mode. In that case an almost 
scalar picture emerges with some polarization effects included to describe the Faraday rotation. In the 
following we will look in more detail what remains of the diffusive picture for Coherent Backscattering 
in a magnetic field. We will not address the role of non-diffusive modes. It is clear that they donot 
influence the angular behavior close to the peak. 

A. Most-Crossed Diagrams 

The most-crossed diagrams associated with Coherent Backscattering can - as usual - be obtained 
by reversing the direction of propagation of the bottom line in the Ladder diagrams. For vector 
waves the bottom polarization indices have to be reversed as well, just like the direction of external 
magnetic fields, in order to remain consistent with the reciprocity relation e (B) = e* (— B) . This 
gives relation (25) and brings us to study the Ladder sum C~ (cu, q, B,0) defined already in Eq. (17). 
Let us first set and q equal to zero and treat them later as perturbative variables. We expect 
that both left and right hand eigenfunction of the eigenvalue of Qq (B) closest to 1 is some linear 
combination of the independent vectors I, <fr and BB (with Gram matrix M). In the notation of 
section IV. A the eigenvalue equations 
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Qo E r * IO = (i-a)E^ 10 ; E^-WCo =(i-a)E 7 .- 

2 2 2 i 

can be rewritten as 



( ri \ 



AM 



, 7i, 7 2 , 7 3 • am-l- =0 



^2 

The bar signifies complex-conjugation. We recall that L~ is the matrix obtained by sandwiching the 
four-rank tensor 1 — Qq between our choice of base vectors. The resulting characteristic equation is, 

det (AM - L~) = . (65) 

We anticipate that A ~ B 2 in leading order of the magnetic field [7]. In this order the solution of 
Eq. (65) is, 

A = i (I" - L" L" ) . (66) 

The second "off-diagonal" term is a subtle vector effect, overlooked in a heuristic theoretical treat- 
ment. The coordinates r,- and of right and left hand eigenfunction become, 

n = 1 h = 1 _ 

r 2 = -L 21 h = ~L 12 (67) 

r 3 = 5 (A — L 31 + L 32 L 2 ^j l 3 = 5 (A — L 13 + L 23 L 12 ^j 

Due to the algebraic complexity it is convenient to express everything in terms of the phase shift a 
(sin a = s) introduced in Eq. (5). Using the explicit form of L~ obtained in the previous section and 
the dimensionless variables /i, ( and w defined in Eq. (7) we obtain , 

A = ^ 2 (2s 2 + l) , (68) 

and, 

r 2 = ifj, ; l 2 = -i[i(2s 2 + 1) 

r 3 = 5A - 2/i 2 ; h = 5A + /i 2 (4s 4 + s 2 - 2) - (w - () sin 2a . 1 ' 

Notice the explicit occurrence of the scattering phase shift a in these formulas, especially the one 
for the the dephasing A. Contrary to what was found earlier for the normal Ladder diagrams £ + , 
the coefficients of the eigenvectors turn out to be complex- valued. For zero phase shift (valid at low 
frequencies) we infer that r,- = 7,-. 

Having found both eigenfunctions and associated eigenvalue for C~ in the diffusive regime, we can 
obtain the algebraic expression for C~ (q, B) in the diffusion approximation. Calculational details 
are left to Appendix D. We obtain the rather simple expression, 

C~ fnORl- 671 |1(-B))(1(B)| 

C (qAB) -TA(5)-(l| K2-(q,ft,B) |r> (70) 
The symmetric form of this expression is - of course - due to the reciprocity relation (23). In 
particular, the minus sign in |1(— B)) is required by the reciprocity principle. 

To obtain the tensor associated with the most-crossed diagrams, polarization indices of the con- 
jugate wave in Eq. (70) have still to be reversed but this can best be done after having added the 
skin layers as well. In the subsections that follow we will find explicitly the diffusion tensor, discuss 
some polarization properties of the most-crossed diagrams, and finally address the shape of the cone 
and the decrease of the enhancement factor. 
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B. Boltzmann Diffusion Tensor in Cone 



Using the eigenfunctions (69) of the four-rank tensor Qq (B) and the matrices S~ and T~ defined 
in Eq. (34), the calculation of the matrix element (1| — 8Q~ (q,0 = 0) |r) in Eq. (70) is possible. To 
extract a diffusion tensor (an intrinsically dynamic transport quantity), a finite frequency difference 
must be incorporated as well. The matrix element for the dynamics is (1| — 8Q~ (0) |r) and 
follows from Eq. (35). The Boltzmann diffusion tensor (B) is defined by bringing the Ladder 
tensor C~ into the form, 

£ "< q ' n ' B) ~-in + q .D 5 (B).q+A(fl)/f (71> 

An anti-symmetric part for (B) has no meaning. We obtain, in terms of the variables (7) and 
the phase-shift a (s = sin a) in Eq. (5), 

(B) = \l (I+</ iso I + 4„iBB) , (72) 



with 



2 / 88 4 907 2 79 \ 1 _ _ . „ 

2 /44 4 97 2 71 \ 1 , ^ . n 
<4ni = H \ —s - —s - — j - -( w - C) sin 2a 

We plotted these corrections as a function of s in Fig. 5. We infer that they clearly deviate from the 
ones of the Boltzmann diffusion constant (B). Two special cases can be discussed. 

The low frequency regime associated with genuine Rayleigh scattering has s ~ (uja) 3 — > 0. In that 
case we find parallel to the field AD^ / D = —6.48 /i 2 and perpendicular to the field AD±/D = 

— 1.75 /i 2 . This implies that diffusion is suppressed in all directions but mainly along the magnetic 
field lines. The last property is in sharp contrast to what was found for the normal incoherent energy 
denoted by C + (as well as for the Beenakker-Senftleben effect) where the suppression was seen to be 
more pronounced perpendicular to B. 

Another special case is when the scatterers are set to resonance: s ~ 1. In that case AD^ / D = 

— 1.6 /i 2 and AD±/D = +0.8 /i 2 . Here, diffusion is even enhanced perpendicular to the field. Due 
to polydispersity this situation may be very difficult to achieve experimentally. Nevertheless, it is a 
nice illustration how the scattering phase shift of the particles may influence macroscopic transport 
phenomena. 

Anisotropy in the diffusion constant may be observed from an anisotropy of the line shape of 
the cone. This will be the topic of the fourth subsection. 



C. Polarization of Cone 

To obtain the most-crossed diagrams at retroflection from the bulk result one must add the skin 
layers as indicated in Eq. (38) and finally transpose the polarization indices of bottom line. We 
will consider only polarization terms linear in the external magnetic field. The normalized lefthand 
eigenfunction then simplifies to, 
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The four-rank tensor in Eq. (70) works out to 

^{liXil+Ixl^Xil+IxliX*!} . 

Next we have to add the skin layers for the incident and emergent light. As has been done for the 
incoherent part the skin layers give rise to a transfer matrix (47), but now with opposite magnetic 
field for the conjugate wave, and a value fx = 1.21 in reflection (Appendix B). Linear in the field it 
reads, 

= A k A k + Ufx [?7iL k A k + T^AkLk] . 

Altogether the polarization matrix for Coherent Backscattering, in the diffusion approximation, 
becomes, 

■ C- (q, B)-MZ k = ^C (B, q) {| A k ) ( A k | + iF [|L k > ( A k | + | A k ) (L k | ]} 

where F = /iRe r/x — fi — 2/is 2 and C (B, q) is the diffusive expression in Eq. (71). In order to get 
the tensor representation for the most- crossed diagrams one must transpose the bottom polarization 
indices k and I and insert q = k + k' (at backscattering this operation leaves the tensor expression 
in the upper formula unchanged). This gives the final result, 

C ijk i (B, k, k') = ^ C (B, k + k') x {A a A jk + iF [L a A kj + A a L kj ]} . (73) 

The scalar factor (2ir/£ 2 )C (B,q) determines the line shape and will be the topic of the next sub- 
section. In particular, we notice that the diffusion approximation decouples polarization and line 
shape. 

With respect to the polarization part it is instructive to consider the experimentally relevant 
situations of linear polarization and circular polarization. Let us first discuss linear polarization 
channels for the case that the magnetic field is perpendicular to the slab (Fig. 7). A polarization 
vector v is incident, w is outgoing. The angle a is defined as the angle between both vectors in the 
direction imposed by the magnetic field using the Lenz Rule. It then follows that, 

w l w k C ljk iv J vi ~ cos 2 a - F sin 2a . (74) 

Obviously the maximum value of this expression is not reached for a = as would happen without 
a magnetic field. The angle with maximum signal is, 

a (max) = arctan 2F ~ —F + O (s 3 ) . (75) 

The rotation angle is linear in the magnetic field and shows up as the combined effect of skin-layer 
and first and last scattering. For Vo > the variable F is strictly negative. In that case the angle 
always obeys the Lenz rule. For low frequencies we find a (max) = 2.21/i radians. At resonance this 
is a (max) = 1.79/i radians. Both the existence of the angle and its sign are in agreement with the 
experiments [11]. 

If the magnetic field is directed along the slab the first polarization modifications due to the 
magnetic field enter in second order of the field and are strictly due to birefringence. Concerning 
circular polarization it is straightforward to show from Eq. (73) that the opposite helicity channel 
has no signal (in diffusion approximation). We found this to be true for orders B 2 in the field as 
well. This agrees with the experiment. 
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D. Dephasing and Line Shape of Cone 



The line shape of Coherent Backscattering is determined by the factor C (q, B) in Eq. (73). In 
an infinite medium it is given by, 

^• B » = A/ f + q .D-(B). q - < 76 » 

Its approximate form for the slab geometry, starting from this bulk result, is obtained from the 
imaging method. For a semi-infinite slab the result was found in Eq. (41), identifying a 2 = X/£ + 

D \\V]\ + ( D l ~ D \\) (B| r qn; 



1 1 + x — exp ( — 2xz n /£) 

R ( q , B) = V 2 1 , (77) 

{ } 2D- x(l + xf K } 



where 



x 2 (B,i,q||) =l 2 X 



X/£ + D^ l + {D- ± -D^ (B r q||) 



, 2 



D 



zz 



determines the dephasing due to both the Faraday effect (A) and a finite backscattering angle de- 
termined by q|| = 2A; sin (0/2) f || ((f). The extrapolation length is now given by Zq ~ 2.1D~ Z and is 
in principle different from the extrapolation length Zq in the Ladder diagrams because the diffusion 
constant is different. The enhancement of the cone relative to the incoherent background can be 
written as, 



1 + x - exp (-2xz^) ^ [ p+ 2zp + 1 polarization cone } 
(2zo + l) x (1 + xf \ D- 2z+ + 1 polarization Ladder J 



Let us discuss these results. First we discuss the suppression of Coherent Backscattering at 
backscattering (8 = 0). Secondly the line shape for finite angles. The suppression at backscattering 
in principle depends on two factors, shown separately in Eq. (78). Both factors depend on both 
direction and magnitude of the magnetic field. If we assume that /i 2 ^> £, w (Appendix C) the 
dimensionless variable ji becomes a one-parameter scaling variable for the enhancement factor. 

The first factor describes the dephasing in multiple scattering due to the Faraday rotation and is 
the main reason why the enhancement factor is suppressed [7]. This contribution dephases similarly 
in angle and in magnetic field, which is also found experimentally up to large fields [11]. 

The second factor in Eq. (78) is due to the different diffusion constants and the different polar- 
ization behavior of Ladder and most-crossed diagrams, giving them different weights. This part is 
only a function of the magnetic field and not of the angle. The polarization will be influenced by 
low orders of scattering as well. At the time of writing no experimental evidence exists for this extra 
factor. We will focus on the first factor. 

An effective medium value for the Verdet constant can be obtained from Eq. (76) by imagin- 
ing the scatterers in a homogeneous medium with Verdet constant V e R (the situation envisaged by 
MacKintosh and John). This defines V e R by, 



X = -V^B 2 £ 2 V eff = (2s 2 + lj ,il(B£) . (79) 
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In particular, using the value for /i at low frequencies given in Appendix C and l/£ = 2fto 4 a 3 
(e — l) 2 / (e + 2) 2 , with / the packing fraction, we get 

Kff ( 3 1 '2e. (80) 



fV \e + 2 

This outcome differs by a subtle factor of a/2 from the effective Verdet constant found from the Dyson 
dispersion law (13). This is due to the fact that the diffusive eigenfunction also changes. For e = 1.15 
Eq. (80) predicts an enhancement over the effective medium value of 1.38. Experimentally a value 
of 1.52 ± 0.15 is obtained [10]. A better agreement between theory and experiment can be obtained 
if a finite value for the phase shift a is adopted since this will augment V e f[. In later experiments [11] 
one introduces the Faraday rotation length £* F > £* in order to distinguish the correlation of circular 
polarization from the correlation in momentum (determined by the usual transport mean free path 
£*). In our theory this length scale does not show up. However, a value £* F /£* > 1 can be translated 
into an enhancement of the effective medium Verdet constant V e f[/ fVo = \J £* F /£* . 

Near resonance (Appendix C) one can relate the variable ji to the "path length" of the wave 
inside the scatterer according to ji = 2VoB L pSLt h so that, 

The Dyson equation gives Vb ys0 n/Vo = —L pa _ t h/£- For strong resonant scattering the path length can 
exceed the mean free path. As a result, the effective Verdet constant will be strongly enhanced by 
resonant scattering. The same mechanism causes the transport speed to go down [23]. A correlation 
between both was suggested by Erbacher, Lenke and Maret [10] but has not been found experimen- 
tally so far. It is well possible that internal spin-flips in Mie scattering destroy this phenomenon 
[38]. 

The line shape of Coherent Backscattering depends on the direction of the magnetic field. When 
the field is perpendicular to the slab, one finds that x 2 ~ A /£D^ + Hence the cone is independent 
of the azimuthal angle <p. In fact this is the standard experimental set-up and no (^-dependence was 
observed. Notice also that in this case the width of the cone is determined by the scattering mean free 
path £, and not the transport mean free path as is known for Mie scattering [28]. This might explain 
why the predicted diminishing of the diffusion constant relevant for the cone due to the magnetic 
field has gone unnoticed so far. 

Experimentally one considered sofar the cone in a magnetic field that is perpendicular to the 
slab. The helicity preserving channel is often preferred because it has the advantage of not giving 
a single-scattering contribution for the incoherent background. Moreover this channel is known to 
have a purely azimuthally isotropic line shape in the absence of a field. It can be verified that 
the polarization impact (73) of the magnetic field on the cone does not change this property. Linear 
polarization channels do not obey this property due to low orders of scattering for which the direction 
of the incident polarization vector is not yet scrambled [39]. 

For arbitrary field direction it follows that 



2 

X ~ 



A /£ + {D^ + (DJ_ - Djj") [l - (B-z) 2 ] cos 2 <p) q 2 
£>,7 + (Dl - D,7) [l - (B-z) 2 



which depends on <p if B is not along the z-axis. In Fig. 8 we display the line shape for the magnetic 
field along the slab. Since we found in Eq. (72) that Dj_ > Dn the cone will more sharpened along 
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the field direction. The azimuthal dependence of the cone with respect to the magnetic field is the 
only experimental way to determine the full diffusion tensor in the most-crossed diagrams. 



VI. FIELD CORRELATION 

In this section we discuss briefly the field correlation function (^Ei (B) Ej (0)^ in transmission. 
The intensity correlation function (/ (B) / (0)) is measurable and has indeed been measured in trans- 
mission. In the Ci-approximation the latter is given by (one plus) the absolute square of the field 
correlation, meaning that a lot of information of the field correlation function can be extracted 
experimentally. 

Using the methods outlined in the previous sections essentially all desired information of the field 
correlation can be obtained analytically. In fact, the way of calculation is quite similar to the one of 
(^Ei (B) Ej (— B)^> which is relevant for Coherent Backscattering. The difference is that the "building 
block" for the Ladder diagrams is not Q~ (B,q) as defined in Eq. (15), but rather 



(82) 



Q s (u, q, B) = J2 G P+f > B ) G * P-§, 0)) • n t (u, B) t* (u, 0) . 

Just like for the cone, the Ladder sum for the field correlation can be written as 
C (u, q, B) = nt (B) t* (0) • {1 + Q s (q) + Q s (q) • Q s (q) + • • •} 
= n.(B,..(0,.E W Gd^l. 

with exactly the same notation as in section V. We restricted again to the two-dimensional subspace 
spanned by left and right eigenvector of Q s (cu, q = 0, B). Let L s (B, q) be the 3x3 matrix containing 
the matrix elements of the four rank tensor 1 — Q s (cu, q, B) with respect to the base {I, BB} . We 
obtain, 

Ln = \ri\ ~ f^3 + L 2 r]2 + 36 - 6 + i?7i6 
L\ 2 = ±i m - 26 

L 21 = ~ 6 

L%2 = 1 + ^! + ^2 ~ i»73 + 6 + |^l6 

£ B 23 = To^i (83) 

£32 = ToiVi ~ \d 

L ls = jotf + 1^2 " f^3 -6 + 6 



L n 2 + h 
10 r n ^ 5' 

L li = h"nl - 1^3 + \it\2 + 6 - |6 + ^6 
£33 = I + jkvl + f^ 2 - I iV3 - |6 + |6 

We will not discuss polarization properties here. The most important result is the remnant of 
long-range diffusion, expressed by the eigenvalue of Q s (q = 0) closest to one. It follows from a 
characteristic equation similar to Eq. (65) that, 

A s = I (L' n - L\ 2 L* 2l ) . (84) 
It can readily be verified that this in fact a complex number. Explicitly, 
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Re A- = i„ 2 (23 2 + l)=^, 

Im A s = --/i 2 sin 2a + -w + -( . (85) 

The real part of the number represents again the dephasing due to the Faraday effect and is a factor 
of 4 less than the dephasing in the crossed diagrams, as was also found in Ref. [10]. The imaginary 
part will give rise to a phase factor in the field-field correlation. As in the expression for Coherent 
Backscattering, 

'e(b)e(o))~ 1 



q D S (B) • q + X s /£ ' 

so that one obtains in transmission from a slab with length L, using formula (41), 

[E(B)E(0)) lV^/£ 



E(0)E(0)) sinh W3A S /^ 



exp(-LV3X s /£) . (86) 



To keep things simple we did not include anisotropy in the diffusion constant (which in fact becomes 
complex- valued as well). Since A s is complex- valued, formula (86) contains an extensive phase factor 
exp (— iQ) with = Im \/3X s L/£ = V e ffBL which has not been considered before. This phase is not 
measurable in the absolute square, although the imaginary part of \/3A s modifies the slope somewhat 
at low fields (see Fig. 8). 

Using the values at low frequencies in Appendix C one finds, 



Re VX 1 
Im a/A^ 



3V6e V B 
2(e-l) 2 (uaf to ' 



1 



4V6e 



5e - 1 V Q B r MB 



2 



U VqB /id 



(87) 



The imaginary part is a factor (loci) 3 smaller than the real part. The low-frequency regime to < a is 
therefore not favorable to observe the phase. Real and imaginary part become comparable when the 
variables /i 2 , ( and w become of the same order of magnitude. This can be achieved when ua 1. 
We emphasize that we deal here with the phase of the complex wave and not with the geometric 
phase associated with the polarization factor. Oscillations in the field correlation (^E x (B) E x (0)^ 
and the intensity correlation (/ (B) / (0)) can originate from the latter due to the fact that the total 
Faraday rotation is always proportional to the total traversed length, even is the path is snakelike and 
not straight [8], quite similar as was found earlier for the Stokes variables in the normal incoherent 
transmission. However, they do not survive in the diffuse regime. 

In Fig. 9 we show field and intensity correlation function. The imaginary part causes a small 
modification in the intensity correlation function at low arguments but does basically not change. 
The field correlation, however, will oscillate due to the presence of the phase 0, and goes through 
zero near a specific value of VqBL. 



VII. CONCLUSIONS 

In this paper we have discussed the impact of magneto-optical effects on multiple elastic scattering 
of light using the framework of transport theory. An inhomogeneous random medium is considered 
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in which the scatterers suffer from the Faraday effect and Cotton-Mouton birefringence and not the 
surrounding medium. The scatterers are modelled by point-like objects, enabling us to get most 
results in closed form. Several aspects of multiple scattering have been addressed. 

We have discussed the complex dispersion law. The complex wavenumber (the real part of which 
gives the effective-medium dielectric constant and the imaginary part the scattering mean free path) 
depends on the direction of propagation with respect to the direction of the magnetic field. Especially 
such dependence on the scattering mean free path has not been discussed elsewhere and causes many 
differences with other work. 

Concerning the average energy density we have calculated the diffusion tensor for the light in 
a magnetic field. We found diffusion to be different along and perpendicular to the field lines. 
The symmetric part follows finally by summing up the Ladder diagrams (equivalent to solving the 
equation of radiative transfer in the diffusion approximation). The anisotropy causes the all-channel 
transmission coefficient of a slab (magneto-conductance) to be different for different orientations of 
the magnetic field. The anti-symmetric part of the diffusion tensor does not feature in the diffusion 
equation but generates a transverse current, perpendicular to magnetic field and energy-density 
gradient. It will be the subject of a separate paper [16]. 

We have considered the Stokes parameters in transmission. For the magnetic field perpendicular 
to the slab we found that a non-zero circular Stokes parameter V persists in diffuse transmission and 
is proportional to the magnetic field strength. The linear Stokes parameters decay exponentially, but 
the presence of the magnetic field causes a rotation in the (£7, Q) plane. If the field is parallel to the 
slab we find that a finite linear polarization remains having the direction of the magnetic field and 
proportional to the square of the field. The circular parameter V now decays exponentially but is 
again accompanied by an oscillatory factor. 

We investigated the most-crossed diagrams, responsible for Coherent Backscattering. The de- 
crease of the factor of two in Coherent Backscattering as a function of the magnetic field (due to the 
modification of the normal reciprocity relations) has been obtained. When the field is parallel to the 
slab we predict the line shape to depend on the azimuthal angle. The diffusion tensor featuring in 
the most-crossed diagrams turned out not to be the same as the one in the Ladder diagrams. Thus 
in a magnetic field one has two different diffusion tensors. In fact, the diffusion tensor for Coherent 
Backscattering was seen to suffer a lot more from the magnetic field and is flattened along the field 
lines. 

The correlation function of the electric field vector with and without external magnetic field has 
been considered theoretically. We obtained a complex-valued diffusive eigenvalue in this case. This 
gives rise to oscillations as a function of the magnetic field. Sofar, only the intensity correlation has 
been measured and not the correlation function of the electric field itself. 

Some questions are very difficult to solve theoretically, and experiments and numerical simulations 
may lead us here. A comparison with the Beenakker-Senftleben effect - as discussed out loud in this 
paper - may help us. For instance, does the diffusion tensor saturate for high magnetic fields? The 
calculation presented in this paper is exact in orders B and B 2 but breaks down when the variable /i 
introduced in this paper becomes comparable to one. Going beyond this regime seems difficult but not 
impossible. Experimentally one has at the moment values maximally equal to ji 1. Concerning 
the dephasing in Coherent Backscattering modifications beyond orders 0(B 2 ) have already been 
reported in multiple scattering experiments [11]. Since modifications in the diffusion tensor have 
been shown to be of the same order we may conclude that this question seems to be resolvable in the 
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near future. In the case of the Beenakker-Senftleben effect the modifications are, in the saturation 
regime - at most a few percents. What happens to the diffusion tensor of light in large fields is 
theoretically not known either. In the Beenakker-Senftleben case the modifications in the saturation 
regime are determined only by the intrinsical asymmetry of the molecules' cross-section. In our case 
the asymmetry of the cross-section is caused by the magnetic field itself. The physical argument 
valid for the Beenakker-Senftleben effect why saturation sets in for large magnetic field, may not 
apply here. 

A qualitative conclusion that one can draw from our calculations is that many properties of 
multiple light scattering in a magnetic field depend explicitly on the scattering phase shift of the 
scatterers. In an experiment one often has a mixture of different kinds of particles and some average 
value should come out. In a monodisperse sample with resonant scatterers however, this notion 
might give rise to rather uncommon features. Usually the phase shift only enters indirectly by means 
of the diffusion constant or the mean free path. 

It is also possible to consider the medium to be magneto-optic rather than the scatterers. We do 
not expect this case to be qualitatively different although no definite conclusions should be drawn for 
this case on the basis of the calculations in the present paper. Even for the aimed situation with the 
magneto-optics inside the scatterers, the model we considered here may be too simplistic. However, 
doing better in the hope of being more close to the experiment will be difficult, if not impossible. 
Moreover, various approximations may overlook the physical phenomena that come out of our model. 
As such we sincerely hope that this paper may serve as a guide for new experiments. 
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APPENDIX A: TENSOR CONVENTIONS 

In general we denote scalars by Roman numbers, second-rank tensors bold, and four-rank tensors 
calligraphic. We will use the convention that a denotes the complex conjugate of a scalar, C* the 
hermitean conjugate and C the complex conjugate of a second-rank tensor C, and the adjoint 
of a four-rank tensor. The second-rank tensor A A denotes the anti-hermitean tensor A — A*. As 
customary, implicit summation over repeated indices is assumed. 

In our notation, the diagram l jT_C_f denotes the four-rank tensor Cijki (in the notation of Ref. [7] 
this would be Cikji)- Two tensors of rank two A and B make up a four-rank tensor (AB) determined 
by (AB) tjkl = A l3 B lk . By definition (AB) C = A • C • B and D (AB) = A* D B* = (A*B*) • D. 
Tensor multiplication occurs as (AB) • (CD) = (A • C) (D • B) . It is easy to show that 

{(AB).(CD)}-F = (AB).{(CD)-F} , 

so that we can forget about the rang of multiplication. A scalar product for second-rank tensors can 
be introduced as 
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(A, B) = Tr A* • B , 

which we will sometimes write as (A | B) . The cyclic property of the trace guarantees that, 

(D,(AB)-C) = (D-(AB),C) = (D| (AB) |C) . 

The adjoint of a four-rank tensor V must obey (A , P • B) = (V^ • A, B^j . It follows that (AB)^ = 
(A*B*) . In general we use the notation, 

<D| C |C) = (D*) kt CtjkiCji , 

where we note that in many physical applications D and C are either hermitean or anti-hermitean. 

APPENDIX B: CALCULATION OF POLARIZATION CONSTANTS 

The constants /,• show up in the calculation of the transfer matrices of the skin layers in reflection 
and transmission due to the Faraday effect and birefringence. They determine an effective length 
with respect to the scattering mean free path I according to 

*-GL« *-l(G)V '-=(£}, • 

1 ' skin \ / skin 



where r = z/ cos 8- m and r' = z' / cos 6 out} both measured from the boundary into the slab. In principle 
they depend on the magnetic field. For our purposes it is sufficient to evaluate these variables for 
zero magnetic field. The constants /,• also depend on the angles 8- m and 6 out . We concentrate on 
waves incident and leaving along the z-axis. 

For Rayleigh scatterers without magnetic field we have 1 = 1*. Let us therefore normalize I = 
£* = 1. For the step length distribution P (z) = exp (—z) - valid far away from the boundaries - it 
follows that f 1 = f 2 = l [36]. The joint distribution of z and z' in reflection is in principle given by 
exp (— z) p (z, z') exp (— z'). Using the diffusion approximation given in Eq. (41) for a semi-infinite 
slab we find, 

P re{ (z,z') = ^— — e~ z [mm (z,z') + z ] e~ z ' , 

z + 1/z 

where the front factor has been chosen such to have normalization 1. Using this distribution we find 
for reflection, (z 0.71) 

^0+1/2 ' " -~o + l/2 ' n r + l/2 

In transmission first (z) and last (V) skin layer are far apart and the joint distribution factorizes into 

-^trans(^)-^trans(^ ) ; 

p / \ _ Z ° Z ~z 

■^trans \^ ) '. T ^ 

zo + l 

In the case of transmission 

/i = ^«1.58; / 2 = ^T^2.16; f 3 = ft . (B3) 

zq + 1 z + 1 
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APPENDIX C: DIMENSIONLESS CONSTANTS 



In this Appendix we relate the dimensionless constants /i, ( and w defined in Eq. (7) and the 
phase shift a in Eq. (5) that occur throughout the paper, to scatterer properties. These variables 
determine the role of magneto-optical effects in multiple scattering. The variable ji signifies the 
modification of reciprocity, determines the Faraday rotation in multiple scattering and is linear in 
the magnetic field; w and ( both determine the birefringence in multiple scattering. The size 
parameter is defined as usual as x = ua with a the radius of the particles. Their volume is given by 
v = 4ira 3 /3. We have, 

" V B 



V B 18e fV B\ 2 9 



u e - 1 )x 3 V u ) 2 e-1 



2MB' 



to 

(CI) 

It is seen that for small size parameters and/or dielectric constants close to one /1 2 ^> (,w, so that 
fi becomes the most important variable. A typical value for ji is (taking the dielectric constant 1.15 
at wavelength 514 nm), 

_ (L73 V B 
^ ~ ~^ 90°/mm/T X 15 T ' 

At very low frequencies i « 1 the variable ji can in fact become quite large. Physically this 
corresponds to a large number of Faraday rotations along one mean free path. The phase shift of 
the scattering matrix at low frequencies is well-known to be 

2e - 1 3 

a = x . 

3e + 2 

The Faraday rotation in the particle is determined by ji. Near a resonance we can find a very 
physical expression. If we neglect the longitudinal field in the scatterer and assume e ^> 1, the 
resonant frequency is Uq = AirT/ev. The path length L p of the light in a particle can be defined as 
the sensitivity of the scattering albedo with respect to absorption [35]. For our model it is inversely 
proportional to the line width UqY and reads L p = 3/2-^/e X 1/cUqT (To have L v ~^> a one must choose 
r 2 < a 2 v^). Then, 

H = 2V BL P . (C2) 

The simplicity of this formula suggests general validity, with only the numerical factor in front 
changing from model to model. The Faraday rotation in the particle would then be proportional 
to the pathlength, independent on any other parameter. For a strong resonance the Faraday effect 
could thus be resonantly enhanced. 

In comparing expression (68) to the one used by Lenke etal. [11] (using a momentum qp 
Knedium-S), leads us to identify fj,y/2 sin 2 a + 1 = qplj\f2. Experimentally one has at most q F £ 2 2 
leading to ji 1. The perturbational treatment of this paper is valid when ji <^1. 



APPENDIX D: DERIVATION OF MOST-CROSSED DIAGRAMS 

In this Appendix we calculate the Ladder diagrams C~ (the minus sign signifying an opposite 
direction of magnetic field for the complex conjugate wave). By the reciprocity principle (25) these 
are related to the most-crossed diagrams. 
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Left and righthand side eigenfunctions relevant for the diffusive regime have been obtained in 
Eq. (69). To avoid subtilities arising from the fact that both eigenfunctions are nearly parallel, we 
prefer the base |r) and |s) = (|r) — |1)) /(ri — l\) and normalize. The matrix N~ is the Gram matrix 
of this new base after normalization. With respect to this base the Ladder diagrams C~ become, for 
finite frequency difference and momentum difference q, 



it(B)t*(-B) • [e-(q,n) + e-( q ,n)-s2-(q,n) + -" 



= B t(B)f(-B). j>> (( N -)- 1 .p- (q ,n). s= -^ y )^(H . 

Notice that we have let the series to start with double scattering in order not to double count single 
scattering. The matrix P~ (q,0,B) is defined by P~ 3 (q,0,B) = (i\ Q~ (q,0,B) \j) and calculation 
yields, 

lq ' ' j ~V° l + A '/ V -i"/«2i + q • D 21 -q -iSl/v 22 + q • D 22 q ) " ^ 
Herein is A' some known constant proportional to B 2 that we shall not need and, 

1 2 
«- (B) = - (l~ 2 + n) = --in (l + s 2 

The entries Vij and D,-j follow straightforwardly from Eqs. (35) and (34). The outlook of Eq. (Dl) 
is quite similar to what was found earlier for N + — P + except that now, for 0, q = 0, the eigenvalue 
with algebraic multiplicity 2 has disappeared in favor of two real-valued eigenvalues A and | + A' 
both with multiplicity 1. The first is proportional to B 2 and the remnant of long-range diffusion. 
From the characteristic equation it is possible to show that - in first order perturbation theory - this 
eigenvalue changes into, 

A(q,0,B) = A - iftun + 3i/c~ft + q- D n -q-3/c~q- D 21 -q=A - (1| 8Q~ (q,0) |r) . (D2) 

The corresponding eigenvector remains approximately the same. In the diffusion approximation we 
ignore all non-diffusive eigenvalues, and we arrive at, 

£-(q,O,B) = nt(B)r(-B).|(|r),|s)).N-.p-(O,0).U.( A(q '° B) ' ^U" 1 ^^} . 

Here the matrix U is the matrix of transformation with the two eigenvectors as column vectors, 

Going through the various matrix multiplications we arrive at the rather simple result that, 

£-( q ,n,B) = „t(B)r(-B). I ^. 

Finally, we can use Eq. (32) to calculate the remaining matrix product. It turns out that, 
nt(B)t*(-B)» = ^ [|I>-Z 1 |*> + Z 2 |BB>] = ^ |1(-B)> . 

This gives the following final outcome for the four-rank tensor C~ (q, B) in the diffusion approxima- 
tion, 

C~ fnORl- 671 Il(-B))(1(B)| 

C (qAB)- TA _ (1| SQ . {qM) |r) (D3) 
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FIGURES 

k(B)/ W 1,(8=0) / l,(B) 
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FIG. 1. Polar plots for the dispersion law at low frequencies, in various directions of propagation in an 
inhomogeneous medium exposed to a magnetic field. The direction of the magnetic field is indicated by 
the arrow. The dashed circle represents the isotropic case, without magnetic field. On the lefthand side 
the wavenumber for the two solutions, right the inverse scattering mean free path. We adopted a packing 
fraction of / = 30%, a dielectric constant e = 1.2 and a magnetic field such that VqB/uj = 0.005. The 
degeneracy of the wavenumber is not resolvable but since the mean free path is very large, the dephasing 
between successive collisions is nevertheless considerable (see Fig. 3). 



k(B)/cj l s (B = 0) / l s (B) 




-10 1 -10 1 



FIG. 2. As in the previous graph but now on resonance (a scattering phase shift of ir/2 has been given 
to the scattering particles). We adopted a size parameter x = 1, dielectric constant e = 1.3, packing fraction 
30% and VqB/u = 0.005. The degeneracy of both wavenumber and scattering mean free path is severe, and 
suffers from Cotton-Mouton birefringence. 
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FIG. 3. The parameter k±£± as a function of orientation with respect to the magnetic field. The same 
parameters have been adopted as in Fig. 1 for Rayleigh scatterers. It can be seen that despite the fact that 
the degeneracy in wavenumber is small, the phase difference (k£) + — (k£)_ that accumulates between two 
succesive collisions is severe, especially for propagation along the magnetic field. 




FIG. 4. Ladder diagrams Cijki and most-crossed diagrams Cijki in a magnetic field. We allow in both 
cases two opposite directions for the magnetic field of the bottom wave. The various possibilities are related 
by reciprocity relations discussed in the text. For pointlike scatterers (whose ^-matrices are indicated by 
crosses) the Ladder diagrams do not depend on p and p', and the Crossed diagrams do not depend on q. 
Bold lines denote the Dyson Green's tensor, dotted lines connect identical particles. The four-rank tensor 
is the "building block" for the Ladder diagrams. 
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FIG. 5. Leading modifications (order B 2 ) to Boltzmann diffusion tensor in a magnetic field. We put 
£ = w = so that the Faraday-rotation parameter fi 2 ~ B 2 determines the modification of the diffusion 
tensors. On the horizontal axis the sinus of the scattering phase shift of one individual scatterer: sin a = 
corresponds to the low-frequency Rayleigh regime; sin a = 1 means on resonance. On top the diffusion 
tensor featuring in the convential Ladder diagrams C + , with equal directions of the magnetic field for wave 
and conjugate wave; On the bottom the diffusion tensor for the Ladder diagrams C~ with opposite direction 
for wave and conjugate wave. This diffusion tensor features in the most-crossed diagrams. Due to the 
absence of reciprocity they are not equal; "para" denotes the diffusion along the magnetic field, "perp" 
perpendicular to the magnetic field. 
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FIG. 6. Anisotropy of the "all channel in all channel out" incoherent transmission (mag- 
neto-conductance) and the forward one-channel incoherent transmission (the incident light is assumed to be 
unpolarized). Both are measured by the relative difference of the Poynting vector for the field lines parallel 
and perpendicular to the slab. Both can be either positive and negative depending on the phase shift a. 




FIG. 7. Slab geometry with the relevant angles to observe the rotation of the polarization vector in 
linear polarization channels of Coherent Backscattering. We find the angle «o with maximum signal to obey 
the Lenz rule for all phase shifts. 
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FIG. 8. Equi-intensity lines for the enhancement factor in Coherent Backscattering in a magnetic field 
for the helicity preserving channel. The magnetic field is directed parallel to the slab, along the horizontal 
axis in this graph. The inset shows the line shapes along (tp = 0, solid) and perpendicular (tp = 7r/2, dashed) 
to the magnetic field. We considered the low-frequency regime (phase shift zero) and took ^ = 0.4. If the 
magnetic field were perpendicular to the slab the equi-intensity lines would be circles with the decrease of 
the enhancement factor being the only impact of the magnetic field. 
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FIG. 9. Intensity (dashed) and field (solid) correlation function. On the horizontal axis the product 

V e ffBL of effective Verdet Constant, magnetic field and slab length. The broad dashed line represents the 

C\ intensity correlation that would have been obtained without imaginary part in A s . For the two other 

1/2 

curves the ratio of imaginary part and real part of A s is 0.8. The real part of the field correlation goes 
through zero as a result of the phase factor. The inset shows the location of the first zero as a function of 
this ratio. For our model this ratio is determined by Eq. (85). 
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